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In Part I we discuss accumulating experimental evidence related to the structure and origin of 
the bosonic spectral function a 2 F(u) in high-temperature superconducting (HTSC) cuprates near 
optimal doping. Some global properties of a 2 F(cu), such as number and positions of peaks, are 
extracted by combining optics, neutron scattering, ARPES and tunnelling measurements. These 
methods give convincing evidence for strong electron-phonon interaction (EPI) with 1 < A < 3 in 
cuprates near optimal doping. Here we clarify how these results are in favor of the Eliashberg-like 
theory for HTSC cuprates near optimal doping. We argue that the neglect of EPI in some previous 
studies of HTSC was based on a number of deceptive prejudices related to the strength of EPI, on 
some physical misconceptions and misleading interpretation of experimental results. 

In Part II we discuss some theoretical ingredients which are necessary to explain the experimental 
results related to pairing mechanism in optimally doped cuprates. These comprise the Migdal- 
Eliashberg theory for EPI in strongly correlated systems which give rise to the forward scattering 
peak. The latter is due to the combined effects of the weakly screened Madelung interaction in 
the ionic-metallic structure of layered cuprates and many body effects of strong correlations. While 
EPI is responsible for the strength of pairing the residual Coulomb interaction (by including spin 
fluctuations) triggers the d-wave pairing. 



PACS numbers: 

I. INTRODUCTION 

In spite of an unprecedentedly intensive experimen- 
tal and theoretical study after the discovery of high- 
temperature superconductivity (HTSC) in cuprates there 
is, even twenty-two years after, no consensus on the pair- 
ing mechanism in these materials. At present there are 
two important experimental facts which are not under 
dispute: (1) the critical temperature T c in cuprates is 
high, with the maximum T c max ~ 160 K; (2) the pair- 
ing in cuprates is d-wave like, i.e. A(k, uo) ~ A s (k,oo) + 
Ad(cj)(cos k x — cosky) with A s < O.IA^. On the con- 
trary there is a dispute concerning the scattering mecha- 
nism which governs normal state properties and pairing 
in cuprates. To this end, we stress that in the HTSC 
cuprates, a number of properties can be satisfactorily 
explained by assuming that the quasi-particle dynam- 
ics is governed by some electron-boson scattering and 
in the superconducting state bosonic quasi-particles are 
gluing electrons in Cooper pairs. Which bosonic quasi- 
particles are dominating in the cuprates is the subject 
which will be discussed in this work. It is known that 
the electron-boson (phonon) scattering is well described 
by the Migdal-Eliashberg theory if the adiabatic param- 
eter Ab = X(ujB/Wb) fulfills the condition Ab <C 1, 
where A is the electron-boson coupling constant, uob is 
the characteristic bosonic energy and Wb is the elec- 
tronic band width. The important characteristic of the 



electron-boson scattering is the Eliashberg spectral func- 
tion a 2 F(k, k', uo) (or its average a 2 F{uo)) which char- 
acterizes scattering of quasi-particle from k to k' by ex- 
changing bosonic energy uo. Therefore, in systems with 
electron-boson scattering the knowledge of this function 
is of crucial importance. There are at least two ap- 
proaches differing in assumed gluing bosons. The first 
one is based on the electron-phonon interaction (EPI) [l[ , 
, , [H , @ where mediating bosons are phonons and 
where the the average spectral function a 2 F{uo) is similar 
to the phonon density of states F p h(uo). Note, a 2 F(uj) is 
not the product of two functions although sometimes one 
defines the function a 2 (uo) = a 2 F (uo) / F (uo) which should 
approximate the energy dependence of the strength of 
the EPI coupling. There are numerous experimental 
evidence in cuprates which support the dominance of 
the EPI scattering mechanism with a rather large cou- 
pling constant 1 < X ep < 3 and which will be dis- 
cussed in detail below. In the EPI approach a 2 F(uo) 
is extracted from tunnelling measurements in conjunc- 
tion with IR optical measurements. We stress again 
that the Migdal-Eliashberg theory is well justified frame- 
work for EPI since in most superconductors the condi- 
tion Ap h < 1 is fulfilled. The HTSC are on the bor- 
derline and it is a natural question - under which con- 
dition can high T c be realized in the non-adiabatic limit 
A p h ~ 1? The second approach @ assumes that EPI is 
too weak to be responsible for high T c in cuprates and it is 
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based on a phenomenological model for spin-fluctuation 
interaction (SFI) as the dominating scattering mecha- 
nism, i.e. it is a non-phononic mechanism. In this ap- 
proach the spectral function is proportional to the imag- 
inary part of the spin susceptibility Imx(k — k',0;), i.e. 
a 2 F(k, k', uo) ~ Imx(k — k',^). NMR spectroscopy and 
magnetic neutron scattering give that in HTSC cuprates 
x(q, uo) is peaked at the antiferromagnetic wave vector 
Q = (it /a, ir /a) and this property is favorable for d-wave 
pairing. The SFI theory roots basically on the strong 
electronic repulsion on Cu atoms, which is usually stud- 
ied by the Hubbard model or its (more popular) deriva- 
tive the t-J model. Regarding the possibility to explain 
high T c solely by strong correlations, as it is reviewed in 
[3], we stress two facts. First, at present there is no vi- 
able theory which can justify these (non-phononic) mech- 
anisms of pairing. Second, the central question in this 
approach is - do models based on the Hubbard Hamilto- 
nian show up superconductivity at sufficiently high criti- 
cal temperatures (T c ~ 100 K) ? A number of numerical 
studies of these models offer a negative answer. For in- 
stance, the sign-free variational Monte Carlo algorithm 
in the 2D repulsive (U > 0) Hubbard model gives no 
evidence for HTSC, neither the BCS- nor Berezinskii- 
Kosterlitz-Thouless (BKT)-like @. At the same time, 
similar calculations show that there is a strong tendency 
to superconductivity in the attractive (U < 0) Hubbard 
model for the same strength of U, i.e. at finite temper- 
ature in the 2D model with U < the BKT supercon- 
ducting transition is favored. Concerning HTSC in the 
t—J model, various numerical calculations such as Monte 
Carlo calculations of the Drude spectral weight [9] and 
high temperature expansion [TO] have shown that there is 
no superconductivity at temperatures characteristic for 
cuprates and if it exists T c must be rather low - few 
Kelvins. These numerical results tell us that the lack of 
high T c (even in 2D BKT phase) in the repulsive (U > 0) 
single-band Hubbard model and in the t — J model is 
not only due to thermodynamical 2D-fluctuations (which 
at finite T suppress and destroy superconducting phase 
coherence in large systems) but it is mostly due to an 
inherent ineffectiveness of strong correlations to produce 
solely high T c in cuprates. These numerical results cer- 
tainly mean that the simple single-band Hubbard, as well 
as its derivative the t-J model, are insufficient to explain 
the pairing mechanism in cuprates and some other in- 
gredients must be included. Having in mind these facts 
there is no room at present for any kind of celebration 
of the victory of non-phononic mechanisms of pairing as 
some prefer to do. 

Since EPI is as a rule strong in oxides, then it is 
plausible that it should be accounted for in cuprates at 
least in the normal metallic state. As it will be argued 
later on, the experimental support for the importance of 
EPI comes from optics, tunnelling, and recent ARPES 
measurements It is worth mentioning that recent 

ARPES activity was a strong impetus for renewed ex- 
perimental and theoretical studies of EPI in cuprates. 



However, in spite of accumulating experimental evidence 
for importance of EPI with X ep > 1, there are occasion- 
ally reports which doubt its importance in cuprates. This 
is the case with recent interpretation of some optical mea- 
surements in terms of SFI only [12[ and with LDA band 
structure calculations pjjj], [Tj], where both claim that 
EPI is negligibly small, i.e. X ep < 0.3. 

The paper is organized as follows. There are two 
parts - in Part I we will mainly discuss experimental 
results in cuprates near optimal doping and minimal 
theoretical explanations which are related to the spec- 
tral function a 2 F(uj) as well to the transport spectral 
function a 2 r F(uo) and how these are related to EPI in 
cuprates. In this work we consider mainly those direct 
one-particle and two-particles probes of low energy quasi- 
particle excitations (by including gap and pseudogap) 
and scattering rates which give informations on the struc- 
ture of the spectral functions a 2 F(k, k', uo) and a 2 r F(uo) 
in systems near optimal doping. These are angle- 
resolved photoemission (ARPES), various arts of tun- 
nelling spectroscopies such as superconductor /insulator/ 
normal metal (SIN) junctions and break junctions, 
scanning-tunnelling microscope spectroscopy (STM) , in- 
frared (IR) and Raman optics, inelastic neutron scatter- 
ing (INS) etc. We shall argue that these direct probes 
give evidence for a rather strong EPI in cuprates as dom- 
inating scattering mechanism of quasi-particles. Some 
other experiments on EPI are also discussed in order to 
complete the arguments for the importance of EPI in 
cuprates. The detailed contents of Part I is the following. 
In Section II we discuss some prejudices related to EPI 
and the Fermi-liquid behavior of HTSC cuprates. We ar- 
gue that any non-phononic mechanism of pairing should 
have very large bare critical temperature T c o ^> T c in 
the presence of the large EPI coupling constant, X ep > 1, 
if its spectral function is weakly momentum dependent, 
i.e. if a 2 F(k,k',cj) « a 2 F(uo). The fact that EPI is large 
in the normal state of cuprates and the condition that it 
must conform with d-wave pairing implies inevitably that 
EPI in cuprates must be strongly momentum dependent. 
In Section III we discuss direct and indirect experimental 
evidence for the importance of EPI and for the weakness 
of SFI in cuprates. These are: 

(A) Magnetic neutron scattering measurements - 
These measurements provide dynamic spin susceptibil- 
ity xCq?^) which is in the SFI phenomenological ap- 
proach @ related to the Eliashberg spectral function, 
i.e. a 2 F s j(k, k', uo) ~ Imx(q = k — k.',uo). We stress 
that such an approach can be qualitatively justified only 
in the weak coupling limit, g s f <C Wb, where Wb is the 
band width. Here we discuss experimental results which 
give evidence for strong rearrangement (with respect to 
uo) of Imx(Q,^) by doping toward the optimal doped 
HTSC [16[. It turns out that in the optimally doped 
cuprates with T c = 92.5 K Imx(Q,^) is drastically sup- 
pressed compared to that in slightly underdoped ones 
with T c = 91 and this is strong evidence for the small- 
ness of the SFI coupling constant. 
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(B) Optical conductivity measurements - From these 
measurements one can extract the transport relaxation 
rate J tr (uj) and indirectly an approximative shape of the 
transport spectral function a 2 F tr (oo). In that respect 
we discuss: (i) the misleading concept concerning the 
relation between the optical relaxation rate J tr (uj) and 
the quasi-particle relaxation rate 7(0;). This (mislead- 
ing) concept has been appearing repeatedly in the last 
twenty years despite the fact that this controversy is re- 
solved many years ago 0, [3, pjj, Q, 0, 0; (ii) 
some methods of extraction of the optical spectral func- 
tion a 2 r F(uj) from optical reflectivity measurements. It 
turns out that the width and the shape of the extracted 
a 2 r F(uj) favor EPI; (iii) the restricted sum-rule for the 
optical weight as a function of temperature which can 
also be explained by strong EPI [20|, [2l|; (iv) good 
agreement with experiments of the temperature depen- 
dence of the resistivity p(T) calculated with the extracted 
a 2 r F(u). Recent femtosecond time-resolved optical spec- 
troscopy in La2- x Sr x Cu04 gives additional evidence for 
importance of EPI and ineffectiveness of SFI [22[ . 

(C) ARPES measurements and EPI - From these 
measurements one can extract the self-energy E(k, uj) 
from which one can extract some properties of 
a 2 F(k, k 7 , uj). Here we discuss the following items: (i) 
appearance of the nodal and anti-nodal kinks in opti- 
mally and slightly under doped cuprates, as well as the 
structure of the ARPES self-energy (X(k, uj)) and its iso- 
tope dependence, which are all due to EPI; (ii) appear- 
ance of different slopes of E(k, uj) at low (uj <C uj p h) and 
high energies (uj ^> uj p h) which can be explained with 
strong EPI; (iii) formation of small polarons in the un- 
doped HTSC which is due to strong EPI - this gives rise 
to phonon side bands which are clearly seen in ARPES 
of undoped HTSC. 

(D) Tunnelling spectroscopy - It is well known that this 
method is of an immense importance in obtaining the 
spectral function a 2 F(uj) from tunnelling conductance. 
In this part we discuss the following items: (i) extraction 
of the Eliashberg spectral function a 2 F(uj) with A = 2 — 4 
from the tunnelling conductance of break-junctions [23[- 
p9^ which gives that the maxima of a 2 F(uj) coincide with 
the maxima in the phonon density of states; (ii) the pres- 
ence of the dip in dl/dV in STM which shows the pro- 
nounced oxygen isotope effect and important role of these 
phonons: (iii) the presence of fine and doping indepen- 
dent structure in I(V) characteristics due to phonon emis- 
sion by the Josephson current in layered HTSC cuprates 
with intrinsic Josephson junctions. 

(E) Phonon neutron scattering measurements - From 
these experiments one can extract the phonon den- 
sity of state F p h(uj) and strengths of the quasi-particle 
coupling with various phonon modes. Here we argue, 
that the large softening and broadening of the half- 
breathing Cu-0 bond-stretching phonon, of apical oxy- 
gen phonons and of oxygen B\ g buckling phonons (in 
LSCO, BISCO,YBCO) cannot be explained by LDA. It 
is curious that the magnitude of softening can be partially 



obtained by LDA but the calculated widths of some im- 
portant modes are an order of magnitude smaller than 
the neutron scattering data show. This remarkable fact 
implies the inadequacy of LDA in strongly correlated sys- 
tems and a more sophisticated many body theory for EPI 
is needed. T his p roblem will be discussed in more details 
in Part II [15]. In Section IV brief summary of the 
Part I is given. Since we are dealing with electron-boson 
scattering in cuprate near optimal doping, then in Sec- 
tion V - Appendix we introduce the reader briefly into 
the Migdal-Eliashberg theory for superconductors (and 
normal metals) where the quasi-particle spectral func- 
tion a 2 F(k, k', uj) and the transport spectral function 
a 2 r F(uj) are defined. 

Finally, at the end of the day one poses a question 
- do the experimental results of the above enumerated 
spectroscopic methods allow a building of a satisfac- 
tory and physically reasonable microscopic theory for ba- 
sic scattering and pairing mechanism in cuprates? The 
posed question is very modest compared with a much 
stringent request for the theory of everything - which 
would be able to explain all properties of HTSC ma- 
terials. Such an ambitious project is not realized even 
in those low-temperature conventional superconductors 
where it is definitely proved that the pairing is due to 
EPI and many properties are well accounted for by the 
Migdal-Eliashberg theory. Let us mention only two ex- 
amples. First, the experimental value for the coherence 
peak in the microwave response cr s (T < T c ) at 17 GHz 
in Nb is much higher than the theoretical value obtained 
by the Migdal-Eliashberg theory [31]. So to say, the 
theory explains the coherence peak at 17 GHz in Nb 
qualitatively but not quantitatively. However, the mea- 
surements at higher frequency ~ 60 GHz are in agree- 
ment with the Migdal-Eliashberg theory [32]. Second, 
the experimental boron (B) isotope effect in MgB^ (with 
T c w 40 K) is much smaller than the theoretical value, 
i.e. a e g P w 0.3 < = 0.5, although the pairing is 
due solely by EPI for boron vibrations [33[. Since the 
theory of everything is impossible in the complex mate- 
rials such as HTSC cuprates in Part I and 77 we shall 
not discuss those phenomena which need much more mi- 
croscopic details and/or more sophisticated many-body 
theory. These are selected by chance: (i) peculiarities of 
the coherence peak in the microwave response cr(T) in 
HTSC cuprates, which is peaked at T much smaller than 
T c , contrary to the case of LTSC where it occurs near 
T c ; (ii) T c dependence on the number of Cu02 in the 
unit cell; (iii) temperature dependence of the Hall coef- 
ficient; (iv) distribution of states in the vortex core, etc. 
However, in a separate paper - Part II [15] we shall dis- 
cuss some minimal theoretical concepts which can explain 
at least qualitatively and semi-quantitatively results re- 
lated to the above enumerated spectroscopic methods. 
Due to the presence of strong correlations and quasi- 
2D electronic structure some of these concepts go be- 
yond the LDA approach. In our opinion at this stage of 
the HTSC physics some important ingredients of the fu- 
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ture theory are already recognized. These are: (1) very 
peculiar quasi-2D ionic-metallic structure with a rather 
weak screening along the c-axis, which is a prerequisite 
for strong EPI; (2) strong Coulomb interaction and cor- 
relations which are responsible for strong magnetism in 
undoped cuprates and for important renormalizations of 
EPI. Since both ingredients belong to the class of strong 
coupling problems, at present there is no quantitative 
theory and therefore we must rely on approximative and 
model theories. Even these approaches allow us qualita- 
tive (and semi-quantitative) explanations of some impor- 
tant properties which are due to the interplay of EPI and 
strong correlations. The latter two cause the appearance 
of momentum dependent EPI - peaked at small transfer 
momenta [2]. Based on such an approach we are able 
to explain (understand) at least qualitatively some very 
puzzling experimental results, for instance: (a) why is 
d-wave pairing realized in the presence of strong EPI? 
(b) why is the transport coupling constant (A tr ) smaller 
than the pairing one A, i.e. X tr < A/3? (c) Why is the 
mean-field (one-body) LDA approach unable to give reli- 
able values for the EPI coupling constant in cuprates and 
how many-body effects help; (d) why is d-wave pairing 
robust in presence of non-magnetic impurities and de- 
fects? (e) why the ARPES nodal and antinodal kinks are 
differently renormalized in the superconducting states? 

In real materials there are numerous experimental ev- 
idence for nanoscale inhomogeneities in HTSC oxides . 
For instance recent STM experiments show rather large 
gap dispersion at least on the surface of BISCO crystals 
|30( giving rise for a pronounced inhomogeneity of the 
superconducting order parameter, i.e. A(k, R) where k 
is the relative momentum of the Cooper pair and R is 
the center of mass of Cooper pairs. One possible reason 
for the inhomogeneity of A(k, R) and disorder on the 
atomic scale can be due to extremely high doping level of 
~ (10 — 20) % in HTSC cuprates which is many orders of 
magnitude larger than in standard semiconductors (10 21 
vs 10 15 carrier concentration). There are some claims 
that high T c is exclusively due to these inhomogeneities 
(of an extrinsic or intrinsic origin) which may increase 
EPI [34]], while other try to explain high T c within the 
inhomogeneous Hubbard or t-J model. In Part II [Hj] we 
argue that the concept of an increase of T c by inhomo- 
geneity is ill-defined, since the increase of T c is defined 
with respect to the average value T c . However, the latter 
quantity is experimentally not well defined and an alleged 
increase of T c by the material inhomogeneity cannot be 
tested at all. 



II. EPI VS NON-PHONONIC MECHANISMS - 
FACTS VS PREJUDICES 

Concerning the high T c in cuprates, one of the central 
questions is - which interaction(s) is(are) responsible for 
strong quasi-particle scattering in the normal state and 
for the superconducting pairing? In the last twenty-two 



years, the scientific community was overwhelmed by all 
kinds of (im)possible proposed pairing mechanisms, most 
of which are hardly verifiable in any material, if at all. 
This trend is still continuing nowadays (although with 
smaller slope), in spite of the fact that the accumulated 
experimental results eliminate all but few. 

A. Fermi vs non-Fermi liquid in cuprates 

After discovery of HTSC in cuprates there was a large 
amount of evidence on strong scattering of quasi-particles 
which contradicts the canonical (popular but narrow) 
definition of the Fermi liquid, thus giving rise to numer- 
ous proposals of the so called non-Fermi liquids, such 
as Luttinger liquid, RVB theory, marginal Fermi liquid, 
etc. In our opinion there is no need for these radical ap- 
proaches in explaining basic physics in cuprates at least 
in optimally, slightly underdoped and overdoped metal- 
lic and superconducting HTSC cuprates. This subject 
will be discussed more in Part II and here we give some 
clarifications related to the dilemma of Fermi vs non- 
Fermi liquid. The definition of the canonical Fermi liquid 
(based on the Landau work) in interacting Fermi systems 
comprises the following properties: (1) there are quasi- 
particles with charge q — ±e, spin s = 1/2 and low-laying 
energy excitations £ k (= ek — ji) which are much larger 
than their inverse life-times, i.e. £ k ^> 1/rk ~ £k/^>- 
Since the level width V = 2/rk of the quasi-particle 
is negligibly small, this means that the excited states 
of the Fermi liquid are placed in one-to-one correspon- 
dence with the excited states of the free Fermi gas; (2) 
at T = K there is an energy level with the Fermi sur- 
face at which £ kF = and the Fermi quasi-particle dis- 
tribution function nj?(^ k ) has finite jump; (3) the num- 
ber of quasi-particles under the Fermi surface is equal to 
the total number of conduction particles (we omit here 
other valence and core electrons) - the Luttinger theo- 
rem; (4) the interaction between quasi-particles are char- 
acterized with a few (Landau) parameters which describe 
low-temperature thermodynamics and transport proper- 
ties. Having this definition in mind one can say that 
if fermionic quasi-particles interact with some bosonic 
excitation (for instance with phonons) and if the cou- 
pling is sufficiently strong, then the former are not de- 
scribed by the canonical Fermi liquid since at energies 
and temperatures of the order of the Debay tempera- 
ture /cbOd(= TujJd) (more precisely ~ 0£>/5), i.e. for 
£k ~ ®d one has r^ 1 > £ k and the quasi-particle picture 
(in the sense of the Landau definition) is broken down. 
In that respect an electron-boson system can be classi- 
fied as a non- canonical Fermi liquid for sufficiently strong 
electron-boson coupling. It is nowadays well known that 
for instance Al, Zn are weak coupling systems since for 
£k ~ ©£> one has r k x <C £ k and they are well described 
by the Landau theory. However, the electron-phonon sys- 
tem is satisfactory described by the Migdal-Eliashberg 
theory and the Boltzmann theory, where thermodynamic 
and transport properties depend on the spectral function 
a 2 F s f(k, k 7 , uo) and its higher momenta. Since in HTSC 
cuprates the electron-boson (phonon) coupling is rather 
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strong and T c is large, i.e. of the order of characteristic 
boson energies (ujb), T c ~ cgb/4, then it is natural that 
in the normal state (at T > T c ) we deal inevitably with 
a strong interacting non-canonical Fermi liquid which is 
at least qualitatively and semi-quantitatively described 
by the Migdal-Eliashberg theory. In order to justify this 
statement we shall in the following elucidate some prop- 
erties in more details by studying optical, ARPES, tun- 
nelling and other experiments. 

B. Prejudice on the limitation of the strength of EPI 
In spite of reach experimental evidence in favor of 
strong EPI in HTSC oxides there was a disproportion in 
the research activity (especially theoretical) in the past, 
since the investigation of the SFI mechanism of pairing 
prevailed in the literature. This retrograde trend was 
partly due to an incorrect statement in [35[ on the pos- 
sible upper limit of T c in the phonon mechanism of pair- 
ing. Since in the past we have discussed this problem 
thoroughly in numerous papers - for the recent one see 
[36( | , we shall outline here the main issue and results only. 

It is well known that in an electron-ion crystal, besides 
the attractive EPI, there is also repulsive Coulomb inter- 
action. In case of an isotropic and homogeneous system 
with weak quasi-particle interaction, the effective poten- 
tial V e ff(\<L,uS) in the leading approximation looks like as 
for two external charges (e) embedded in the medium 
with the total longitudinal dielectric function £tot(k, uj) 
(k is the momentum and uj is the frequency) [37[, i.e. 



V eff (k,uj) = 



Aire 2 



e to t(Kuj) k 2 e tot (\z,uj)' 



(1) 



In case of strong interaction between quasi-particles, the 
state of embedded quasi-particles changes significantly 
due to interaction with other quasi-particles, giving rise 
to V e ff(\a,uj) 7^ Aire 2 /k 2 s t ot(k, uj). In that case V e ff 
depends on other (than e tot (k, a;)) response functions. 
However, in the case when Eq.(pQ) holds, i. e. when 
the weak-coupling li mit is realized, T c is given by T c w 
cjexp(-l/(A ep -/i*) [13, [13). Here, X ep is the EPI cou- 
pling constant, uj is an average phonon frequency and /jl* 
is the Coulomb pseudo-potential, /x* = /x/(l + /iln Ep/O) 
(Ef is the Fermi energy). The couplings X ep and \i are 
expressed by £t t(k, uj = 0) 

M -A ej) = (AT(0)T/ e// (k, W = 0)) 



N(0) 



JO 



kdk 



Aire 2 



2k 2 F k 2 e tot (k,uj = 0)' 



(2) 



where N(0) is the density of states at the Fermi sur- 
face and kp is the Fermi momentum - see more in [1]. 
In [35[ it was claimed that lattice stability of the sys- 
tem with respect to the charge density wave forma- 
tion implies the condition £t t(k, uj = 0) > 1 for all 
k. If this was correct then from Eq.([2]) it follows that 
ji > A ep , which limits the maximal value of T c to the 
value T c max « £ F exp(-4 - 3/X ep ). In typical metals 



Ef < (1 — 10) eV and if one accepts the statement in 
jaai, i.e. that X ep < fi < 0.5, one obtains T c ~ (1 - 10) 
K. The latter result, if it would be true, means that EPI 
is ineffective in producing not only high-T c superconduc- 
tivity but also low-temperature superconductivity (LTS). 
However, this result is apparently in conflict first of all 
with experimental results in LTSC, where in numerous 
systems \i < X ep and X ep > 1. For instance, X ep ^2.6 
is realized in PbBi alloy which is definitely much higher 
than /i(< 1). 

Moreover, the basic theory tells us that e to t(k 7^ 0, w) is 
not the response function [37[ . Namely, if a small external 
potential 5V ex t(k, uj) is applied to the system it induces 
screening by charges of the medium and the total poten- 
tial is given by 5V to t(K^) = 5V ex t{K uj)/stot(K u) which 
means that 1/etot (k, uj) is the response function. The lat- 
ter obeys the Kramers-Kronig dispersion relation which 
implies the following stability condition: l/s to t(k,uJ = 
0) < 1 for k ^ 0, i.e. either e tot (k + 0,uj = 0) > 1 
or £tot(k 7^ 0,o; = 0) < 0. This important theorem in- 
validates the above restriction on the maximal value of 
T c in the EPI mechanism. We stress that the condition 
£tot(k ^ 0, uj = 0) < is not in conflict with the lat- 
tice instability. For instance, in inhomogeneous systems 
such as crystal, the total longitudinal dielectric function 
is matrix in the space of reciprocal lattice vectors (Q), 
i.e. £tot(k + Q, k + Q', a;), and £tot(k, uj) is defined by 
£^J(k, uj) = £^,t(k + 0, k + 0, uj). It is well known that 
in dense metallic systems with one ion per cell (such 
as metallic hydrogen ) and with the electronic dielectric 
function £ e z(k, 0), one has [38^ 



£tot (k,0) 



£ e z(k,0) 



l-l/e cZ (k,0)G ep (k) 



(3) 



At the same time the frequency of the longitudinal 
phonon uJi(k) is given by 



£e/(k,0) 



[l-eei(k,0)G ei ,(k)], 



(4) 



where Q 2 is the ionic plasma frequency, G ep is the local 
(electric) field correction - see Ref. [38| . The real condi- 
tion for lattice stability requires that uJ 2 (k) > 0, which 
implies that for e e z(k, 0) > one has e e z(k, 0)G ep (k) < 1. 
The latter condition gives automatically e tot (k, 0) < 0. 
Furthermore, the calculations [38[ show that in the metal- 
lic hydrogen crystal, Etotik-, 0) < for all k 7^ 0. More- 
over, the analyzes of crystals with more ions per unit 
cell |38| gives that s to t(k^0,0) < is more a rule 
than an exception - see Fig. [TJ The physical reason for 
£tot(k 7^ 0,0) < are local field effects described above 
by G ep (k). Whenever the local electric field Ej oc act- 
ing on electrons (and ions) is different from the aver- 
age electric field E, i.e. Ej oc 7^ E, there are correc- 
tions to Etot (k, 0) (and to £ e (k, 0)) which may lead to 
£ £o£ (k,0)<0. 



6 



-! 



-2 



6"'(P,0) >-V.AI 

^S&zJ P 




\Pb 



FIG. 1: Inverse total static dielectric function e~ 1 (p) for nor- 
mal metals (K, Al, Pb and metallic H) in p = (1, 0, 0) direc- 
tion. G is the reciprocal lattice vector. 



C. Is a non-phononic pairing realized in HTSC? 

Regarding EPI one can pose a question - whether it 
contributes significantly to d-wave pairing in cuprates? 
Surprisingly, despite numerous experiments in favor of 
EPI, a number of researchers still believe that EPI is ir- 
relevant for pairing [6[. This belief is mainly based: (i) 
on the, previously discussed, incorrect lattice stability 
criterion, which implies small EPI; (ii) on the well estab- 
lished experimental fact that d-wave pairing is realized 
in cuprates [4l[, which is believed to be incompatible 
with EPF Having in mind that EPI in HTSC is strong 
with 1 < X ep < 3 (see below), we assume for the mo- 
ment that the leading pairing mechanism in cuprates, 
which gives d-wave pairing, is due to some non-phononic 
mechanism, like the exitonic one, with the high energy 
gluing boson (Q np h ^> Wph) an d with the bare critical 
temperature T c o and look for the effect of EPF If EPI 
is isotropic, like in most FTSC materials, then it would 
be very detrimental for d-wave pairing - the pair break- 
ing effect. In the case of dominating isotropic EPI in 
the normal state and the exitonic-like pairing, then near 
T c the linearized Eliashberg equations have an approxi- 
mative form for weak non-phonon interaction (with the 
characteristic frequency Q n ph) 



The above analysis tells us that in real crystals 
£tot(k, 0) can be negative in the large portion of the Bril- 
louin zone giving rise to X ep — n > in Eq.([2]). This 
means that the dielectric function e to t does not limit T c 
in the phonon mechanism of pairing. This result does 
not mean that there is no limit on T c at all - see more 
in [36[ and references therein. We mention in advance 
that the local field effects play important role in HTSC 
oxides, due to their layered structure with very unusual 
ionic-metallic binding, thus giving rise to large EPI - 
see more in the subsequent sections. It is pertinent to 
note that one of the author of [35[ recognizes the pos- 
sibility £tot(k, 0) < and in [39[ even makes interesting 
proposals for compounds with large EPI and T c > 100 
K, while the other author [40] still ignores rigors of scien- 
tific arguments and negates importance of EPI in HTSC 
cuprates. 

In conclusion we point out that there are no theoretical 
and experimental arguments for ignoring EPI in HTSC 
cuprates. To this end it is necessary to answer several im- 
portant questions which are related to experimental find- 
ings in HTSC cuprates (oxides): (1) if EPI is responsible 
for pairing in HTSC cuprates and if superconductivity 
is of d — wave type, how are these two facts compati- 
ble? (2) why is the transport EPI coupling constant X tr 
(entering resistivity) much smaller than the pairing EPI 
coupling constant X ep (> 1) (entering T c ), i.e. why one 
has A* r (^ 0.4 - 1.2) <C A ep (- 2 - 4)? (3) is high T c 
possible for a moderate EPI coupling constant, let say 
for X ep < 1, and under which conditions? (4) if EPI is 
ineffective for pairing in HTSC oxides, inspite of X ep > 1, 
why it is so? 



Z(w n )A(k,w n ) =nT c J2 ^K P fc(k,q,n,n0-f^p 



(5) 



For pure d-wave pairing one has V np h (k, q, n, ra ) = 
V nph Q(Q nph - \uj n \)Q(Q nph - \u) n >\) x (cosk x - 
cosk y )(cosq x - cosq y ) and A(k,cj n ) = A d Q(Q nph - 
\u) n \)(cos k x — cos k y ) which gives the equation for T c - 
see [l| 



In- 



t c 



1 



- ep 



2nT r 



)• 



(6) 



Here ^ is the di-gamma function. At temperatures near 
T c one has T ep w 2i\X ep T c and the solution of Eq. (|6]) is 
approximately T c w T c o exp{— A ep }, which means that for 
T c max - 160 K and X ep > 1 the bare T c0 due to the non- 
phononic interaction must be very large, i.e. T c o > 500 
K. 

Concerning other non-phononic mechanisms, such as 
the SFI one, the effect of the isotropic EPI in the frame- 
work of Eliashberg equations was studied numerically in 
[42| . The latter is based on Eas. (l42M4p in Appendix A. 
with kernels 



^kp(^n) = Akp(^n) + Ajfp(^n) 



AkpOn) = A^(^n) - Ak^n), 



(7) 



(8) 



where X^ p (iv n ) is taken in the FFEX approximation [43j]. 
The calculations [42[ confirm the very detrimental effect 
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of the isotropic EPI on the d-wave pairing due to SFI. For 
the bare SFI T c0 ~ 100 K and X ep > 1 the calculations 
give very small T c <C 100 K. These results tell us that 
a more realistic pairing interaction must be operative in 
cuprates and that EPI is strongly momentum dependent 
[451 ]. Only in that case is strong EPI conform with d- 
wave pairing, either as its main cause or as a supporter 
of a non-phononic mechanism. In Part II we shall argue 
that the strongly momentum dependent EPI is the main 
player in cuprates providing the strength of the pairing 
mechanism, while the residual Coulomb interaction and 
SF, although weaker, trigger it to d-wave pairing. 



III. EXPERIMENTAL EVIDENCE FOR 
STRONG EPI AND WEAK SFI 

In the following we discuss some important experi- 
ments which give evidence for strong EPI in cuprates. 
Before doing it; we shall discuss some magnetic neutron 
scattering measurements in cuprates whose results are 
against the SFI mechanism of pairing. The experimen- 
tal results related to the pronounced imaginary part of 
the susceptibility Jmx(k, k z , uj) at the AF wave vector 
k = Q = (7r, 7r) were interpreted in a number of papers 
as a support for the SFI mechanism for pairing [6]. We 
briefly explain inadequacy of such an interpretation. 



A. Magnetic neutron scattering and the spin 
fluctuation spectral function 

A. SFI affects T c very little 

Before discussing experimental results in cuprates on 
the imaginary part of the spin susceptibility 7mx(k, uj) 
we point out that in theories based on spin fluctuations 
the effective pairing potential V^/(k, a;), which is repul- 
sive, is assumed in the form [6] 



F s/ (q,^ + i0 + ) =g 2 sf 



dv /mx(q, v + i0 + ) 



v — UJ 



(9) 



This form of V s f can be theoretically justified in the weak 
coupling limit (U <C Wb) only. This mechanism of pair- 
ing could be effective in cuprates only if the spin suscep- 
tibility (spectral function) Imx(q, uj) is strongly peaked 
at the AF wave vector Q = (7r/a,7r/a). What is the ex- 
perimental situation? The breakthrough came from mag- 
netic neutron scattering experiments on Y Ba2C r U2 > 0§+ x 
by Bourges group (161 ] . They showed that Imx^(q, uj) 
(the odd part of the spin susceptibility in the bilayer sys- 
tem) at q = Q = (7r, 7r) is strongly dependent on the hole 
doping as it is shown in Fig. O 

By varying doping there is a huge rearrangement of 
Imx^^iQiUj) in the frequency interval which is impor- 
tant for superconducting pairing, let say 5 meV < uj < 60 
meV as it is seen in the last two curve in Fig. [2jtop). 
At the same time there is only a small variation of the 



Normal State, 100 K, Q=(m} 




c ISO 

3 

t m 

i 
< 

! « 



. 1 1 YBa CuO 

m ; Hi 



or 



. L . .... •• •• 



n » N « St « 
Eiwgy(iiteV) 

SC State, 5K, Q={a) 




z « 

fc m 
I « 

< . 

! 

I 

m 

*S5 



J. 



/ \ T »St K 



VBaCuO k 



YBiCtiO 



fit t 



T £ = 52.SK 



T c «HK 



-Ijt 



it: 

■ 
t 



14 21 » 44 SO H 
Ewgf (<ntV) 



FIG. 2: Magnetic spectral function Imx^ (k, cj) in 
y^Ba^Ci^Oe+x- (Top) In the normal state at T = 100 K and 
at Q = (tt, tt). 100 counts in the vertical scale corresponds to 
Xmax ~ 350/i^/eV. (Bottom) In the superconducting state 
at T = 5 K and at Q = (tt, tt). From Ref. 0. 
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corresponding critical temperature T c \ For instance, in 
the underdoped YBa^Cu^O^^ crystal Imx^ od ^(Q, u), 
and 5(Q) = N(fi)g 2 sf du Im X {odd) (Q^) ~ A s/ • (u), 
is much larger than that in the near optimally doped 
YBCL2CU3O6.97, i.e. 6 , 6 .92(Q) > 5 , 6.97(Q), although the 
difference in the corresponding critical temperatures T c 
is very small, i.e. J^ 6 ' 92 ^ = 91 K (in YBa2CuzO§^2) and 
j,j6.97) _ g2 ^ ^ n YBa2CusOe.gr). This pronounced 
rearrangement and decrease of Imx^ odd ^(Qj ^ by dop- 
ing, but a negligible change in T c in YBCO is clearly 
seen in Fig. [2ttop), which is strong evidence against the 
SFI mechanism of pairing. These results in fact mean 
that the SFI coupling constant A s *^(~ g 2 ^) is small, i.e. 

Aj xp) <1, and the SFI pairing mechanism is ineffective 
in cuprates. We stress that in the phenomenological the- 
ory of the SFI pairing [6[ , an unrealistically large coupling 
o s f > 0.7 eV was assumed which gives A s/ - 2. The lat- 
ter value cannot be justified neither experimentally nor 
theoretically. Let us add that the anti- correlation be- 
tween the decrease of Imx(Q,o;) and increase of T c by 
increasing doping toward the optimal value is also present 
in the NMR spectral function Iq = lim^^o I m x(Q> w)/u) 
which determines the longitudinal relaxation rate 1/Ti - 
see 0. This result additionally disfavors the SFI model 
of pairing [6], i.e. the strength of pairing interaction is 
little affected by SFI. As we shall discuss below the role 
of SFI together with the stronger direct Coulomb inter- 
action is to trigger d-wave pairing. 

A less direct argument for smallness of the SFI cou- 
pling constant, i.e. g s f < 0.2 eV and X s ^ ~ 0.2 — 0.3 
comes from other experiments related to the magnetic 
resonance peak in the superconducting state, and this 
will be discussed next. 

B. Ineffectiveness of the magnetic resonance peak 

In the superconducting state of optimally doped 
YBCO and BISCO, Imx(Q,o;) is significantly sup- 
pressed at low frequencies except near the resonance en- 
ergy oj res « 41 meV where a pronounced narrow peak 
appears - the magnetic resonance peak. We stress that 
there is no magnetic resonance peak in LSCO sand con- 
sequently one can question the importance of the reso- 
nance peak in the scattering processes. The relative in- 
tensity of this peak (compared to the total one) is small, 
i.e. Jo ~ (1 — 5)% - see Fig [2] (bottom). In underdoped 
cuprates this peak is present also in the normal state as it 
is seen in Fig[2](top). After the discovery of the resonance 
peak there were attempts to relate it: (i) to the origin of 
the superconducting condensation energy and (ii) to the 
kink in the energy dispersion or the peak-dimp structure 
in the ARPES spectral function. In order that the prop- 
erty (i) holds it is necessary that the peak intensity io is 
small [46]. Io is obtained by equating the condensation 
energy E con with the change of the magnetic energy in 
the superconducting state, i.e. 5E mag « 4:IoE mag , where 

E con « A^(0)A 2 /2 



Fmag = J J J i 1 ~ COS k x ~ cos k y )S(k, u) . (10) 

By taking 2 A w 4T C and the realistic value N(0) ~ 
1/(10 J) ~ 1 states jeV • spin, one obtains Io ~ 
10 _1 (T C /J) 2 ~ 10 -3 . However, such a small intensity 
cannot be responsible for the anomalies in ARPES and 
optical spectra since it gives rise to small coupling con- 
stant X res for the interaction of holes with the resonance 
peak, i.e. X res « (2I N(0)g 2 f /u; res ) <C 1. Such a small 
coupling does not affect superconductivity. Moreover, by 
studying the width of the resonance peak one can extract 
the SFI coupling constant g s f. Thus, the magnetic res- 
onance disappears in the normal state of the optimally 
doped YBCO, which can be qualitatively understood by 
assuming that its broadening scales with the resonance 
energy uo resi i.e. 7 res < uo res , where the line- width is 
given by 7 res = 47r(N(0)g sf ) 2 uj res 0. This limits to 
g s f < 0.2 eV. We stress that the obtained g s f is much 
smaller than the one assumed in the phenomenological 
spin-fluctuation theory Q where g s f ~ 0.6 — 0.7 eV, 
but much larger than in [46] (where g s f < 0.02 eV). 
The smallness of g s f comes out also from the analysis 
of the antiferromagnetic state in underdoped metals of 
LSCO and YBCO [47[, where the small magnetic mo- 
ment /i(< 0.1 n B ) points to an itinerant antiferromag- 
netism with small coupling constant g s f < 0.2 eV. The 
conclusion is that the magnetic resonance in the opti- 
mally doped YBCO is a consequence of the onset of su- 
perconductivity and not its cause. 

There is also a principal reason against the pairing 
due to the resonance peak at least in optimally doped 
cuprates. Since its intensity near T c is vanishingly small, 
though not affecting pairing at the second order phase 
transition at T c , then if it would be the origin for su- 
perconductivity the phase transition at T c would be first 
order, contrary to experiments. Recent ARPES experi- 
ments give evidence that the magnetic resonance cannot 
be related to the kinks in ARPES spectra [jj], [49^ - see 
the discussion below. 

We shall argue below that despite its smallness, spin 
fluctuations can, together with other contributions of 
the residual Coulomb interaction, trigger d-wave pair- 
ing, while the strength of pairing is due to EPI which is 
peaked at small transfer momenta - see more below and 
in 0, 0. 



B. Optical conductivity and EPI 

Optical spectroscopy gives information on optical con- 
ductivity cf{uo) and on two-particle excitations, from 
which one can indirectly extract the transport spec- 
tral function a 2 r F(uo). Since this method probes bulk 
sample (on the skin depth), contrary to ARPES and 
tunnelling methods which probe tiny regions (10 — 15 
A) near the sample surface, this method is very in- 
dispensable. However, a(u) is not a directly measured 
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quantity but it is derived from the reflectivity R(ou) — 

2 

i) 



with the transversal di- 



(y/eu(u) - l)/(^£u(ou) 

electric tensor £u(ou) = £„ >00 + £u,iatt + Airier u{ou) /ou. 
Here, £u i0 c is the high frequency dielectric function, 
£u,iatt describes the contribution of the lattice vibrations 
and (Ju(uo) describes the optical (dynamical) conductiv- 
ity of conduction carriers. R{ou) was usually measured in 
the limited frequency interval o; m i n < ou < ou max . There- 
fore, some physical modelling for R(ou) is needed in order 
to guess it outside this range - see more in reviews 

0. This was the reason for numerous inadequate inter- 
pretations of optic measurements in cuprates, as well as 
the misconceptions and misinterpretations that will be 
uncover below. An illustrative example for this claim is 
large dispersion in the reported value of ou p i - from 0.06 
to 25 eV, i.e. almost three orders of magnitude - see dis- 
cussion in [Hcj. This tells us also that in some periods 
science suffers from a lack of rigorousness and objective- 
ness. However, it turns out that IR measurements of 
R{ou) in conjunction with elipsometric measurements of 
£u(ou) at high frequencies allows reliable determination 
of (j(ou). 

1. Transport and quasiparticle relaxation rates 
The widespread misconception in studying the quasi- 
particle scattering in cuprates was an ad hoc assumption 
that the transport relaxation rate 7^ r (cc;) is equal to the 
quasi-particle relaxation rate 7(0;), in spite of the well 
known fact that j tr (w) 7^ 7(^0 [ill- This incorrect as- 
sumption led to the abandoning of EPI as relevant scat- 
tering mechanism in cuprates. Although we have dis- 
cussed this problem several times before, we want to do 
it again, since the correct understanding of the scatter- 
ing mechanism in cuprates will take us forward in under- 
standing of the pairing mechanism. 

The dynamical conductivity cr(ou) consists of two parts, 

1. e. a(ou) = a inter (ou) + a intra (ou) where a inter (ou) de- 
scribes interband transitions which contribute at higher 
frequencies, while a intra (ou) is due to intraband transi- 
tions which are relevant at low frequencies ou < 1 eV. (In 
IR measurements the frequency is usually given in cm -1 , 
where the following conversion holds: 1cm -1 = 29.98 
GHz = 0.123985 meV = 1.44 K.) The experimental 
data for cr(ou) = o~\ + i&2 in cuprates are usually pro- 
cessed by the generalized (extended) Drude formula [17[ , 

& EH, M 



(j(ou) 



1 



1 



47r j tr (ou) - ioum tr {ou) I 'rrioo ou tr (ou) ' 



(11) 



which is a useful representation for systems with sin- 
gle band electron-boson scattering which is justified in 
HTSC cuprates - see the discussion below. (The useful- 
ness of introducing the optic relaxation ou tr (u) will be 
discussed in Appendix B.) Here, i = a,b enumerates the 
plane axis, ou pi j tr (uj,T) and m op (ou) are the electronic 
plasma frequency, the transport (optical) scattering rate 
and the optical mass, respectively. Very frequently, 
the quantity j* r (ou : T) = % r (cj, T^m^ / m tr {ou)) = 



Im a(ou)/ou Re a(ou) [5l|, which is determined from the 
half- width of the Drude-like expression for cr(cj), was an- 
alyzed since it is independent of ou p . In the weak coupling 
limit X ep < 1, the formula for conductivity given in Eqs. 
(j64l!67[ ) can be written in the form of Eq. ([TT]) where j tr 
reads]lH-[l9l 



PC 

^ L 

j Jo 



di/a? P ,F,(i/)[2(l 



■2n B (i/)) 



V to + V , . to — V . 
-2 n B {u + v) + n B {u-v)\. (12) 

OU OU ou 

Here ub{ou) is the Bose distribution function, 
explicit form of the transport mass m tr (ou) see 
03, HL 00 l n the presence of impurity scattering one 
should add jl™ p to j tr . It turns out that Eq. ([T2|) holds 
within a few percents also for large X ep (> 1). Note, that 
a tr 7^ afFi(v) and the index I enumerates all scat- 

tering bosons - phonons, spin fluctuations, etc. For com- 
parison, we give the quasi-particle scattering rate 7(0;, T) 



00 

,T) = 2tt J dva 2 F(y){2n B (y) 



70 



+ n F {y + ou) + n F {y - ou)} + Y 



(13) 



where n F is the Fermi distribution function. By com- 
paring Eq.([T3]) and Eq. ([T2]h it is seen that j tr and 7 are 
different quantities, j tr 7^ 7, i.e. the former describes the 
relaxation of Bose particles (electron-hole pairs) while the 
latter one the relaxation of Fermi particles. This differ- 
ence persists also at T = K where one has (due to 
simplicity we omit in the following summation over I ) 

da 



2tt 



/ <&/(. 
Jo 



to - v)a 2 tr (v)F(v) (14) 



and 



j(uj) = 2tt / dvo?{y)F(y). 
Jo 



(15) 



In the case of EPI, the above equations give that 
j ep (ou) = const for ou > ou™* while 7^(0;) (as well as ^ r ) 
is monotonic growing for ou > ou™*, where ou™* is the 
maximal phonon frequency. This is clearly seen by com- 
paring 7(0;, T), 7 £r (cj, T) and 7^ which are calculated for 
the EPI spectral function al p (ou)F p h(ou) extracted from 
tunnelling experiments in YBCO (with ou™* ~ 80 meV) 
H3 - see Fig. [S 

The results shown in Fig. [3] clearly demonstrate the 
physical difference between two scattering rates j ep and 
7^. It is also seen that 7^(0;, T) is more a linear function 
of ou than j tr (ou,T). From these calculations one con- 
cludes that the quasi-linearity of j tr (ou,T) (and 7^) is 
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FIG. 3: (a) Scattering rates 7(0;, T), 7 tr (o;, T) and 7^ - from 
top to bottom, for the Eliashberg function in (b). From [l8[. 
(b)Eliashberg spectral function a^ r) (uj)F v h(uo) obtained from 
tunnelling experiments on break junctions [23]. Inset shows 
7^ r with (full line) and without (dashed line) interband tran- 
sitions protect [l[. 



not in contradiction with the EPI scattering mechanism 
but it is in fact a natural consequence of EPI. We stress 
that such behavior of j ep and 7^, shown in Fig. [3l is in 
fact not exceptional for HTSC cuprates but it is generic 
for many metallic systems, for instance 3D metallic ox- 
ides, low temperature superconductors such as A/, Pb, 
etc. - see more in [l|, Q. 

Let us discuss briefly the experimental results for R(uj) 
and 7^ r (cc;,T) and compare these with theoretical pre- 
dictions obtained by using a single band model and 
al p (co)Fph(co) from tunnelling data with the EPI cou- 
pling A = 2 [23]. In the case of YBCO the agreement 
between measured and calculated R(uo) is very good up 
to frequencies uo < 6000 cm -1 which confirms the im- 
portance of EPI in scattering processes. For higher fre- 
quencies, where a mead infrared peak appears, it is nec- 



essary to account for interband transitions [l|. In opti- 
mally doped Bi2Sr2CaCu20e [52j] the experimental re- 
sults for 7^(0;, T) are explained theoretically by assum- 
ing that the EPI spectral function a 2 ep (uo)F(uo) ~ F p h(uo), 
where F p h(uo) is the phononic DOS in BISCO while 
^epi 00 ) ~ u 1 ' 6 ! ^ = 1-9 an d lim ~ 320 era -1 - see Fig. 
Ufa). The agreement is rather good. At the same time 
the fit of jt r (uo,T) by the marginal Fermi liquid fails as 
it is evident in Fig. 2Kb). 

Now we will comment on the so called pronounced lin- 
ear behavior of j tr (w, T) (and 7^ r (c<;, T)) which served in 
the past for numerous inadequate conclusions. We stress 
that the measured quantity is reflectivity R(uo) and de- 
rived ones are cr(u)), J tr (w, T) and m tr (oo), which are very 
sensitive to the value of the dielectric constant Sqq . This 
is clearly demonstrated in Fig. [5] for Bi2212 where it is 
seen that 7 tr (o;,T) (and 7^(0;, T)) for Soo = 1 is linear 
up to much higher uo than in the case ^ !• 

In some experiments [53|, [HH ^ tr (uo, T) (and 7* r (k;, T)) 
is linear up to very high uo which means that the ion 
background and interband transitions (contained in e^) 
are not properly taken into account since it is assumed 
too small Soq. The recent elipsometric measurements on 
YBCO [55[ give the reliable value for w 4 — 6. The 
latter gives rise to a much less spectacular linearity in the 
relaxation rates than it was the case immediately after 
the discovery of HTSC cuprates. 

Furthermore, we would like to comment on two points 
concerning cr, 7 £r , 7 and their interrelations. First, the 
parametrization of cr(uo) with the generalized Drude for- 
mula in Eq. (fTT]) and its relation to the transport scat- 
tering rate j tr (uo, T) and the transport mass m tr (uo, T) is 
useful if we deal with electron-boson scattering in a single 
band problem. In [19] it is shown that <j(uo) of a two-band 
model with only elastic impurity scattering can be repre- 
sented by the generalized (extended) Drude formula with 
uo and T dependence of effective parameters 7^/^(0;, T), 
rrifrf(u,T) despite the fact that the inelastic electron- 
boson scattering is absent. To this end we stress that 
the single-band approach is fully justified for a number 
of HTSC cuprates such as LSCO, BISCO etc. Second, 
at the beginning we said that j tr (uo,T) and j(uo : T) are 
physically different quantities and it holds j tr (uo,T) ^ 
7(0;, T). In order to give the physical picture and quali- 
tative explanation we assume that a 2 r F(v) ~ a 2 F(v). In 
that case the renormalized frequencies, the quasi-particle 
one uo(uo) = Z(uo)uo = uo — Yi(uo) (for the definition of Z{uo) 
see Appendix A.) and the transport one ub tr (u) - defined 
above, are related and at T = 0, they are given by pj}, 

M 



1 r u 

uo tr (uo) = - / duj'2Q((jj'). (16) 
^ Jo 

It gives the relation between 7 £r (c^) and 7(0;), m tr (uo) 
and m*(uj) respectively 

7trH = - j" dLj'-iiw') (17) 

w Jo 
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FIG. 4: (Top) Experimental transport scattering rate 7^ r 
(solid lines) for BISCO and the theoretical curve by using 
Eq. ([64]) and transport mass m* tr with a 2 F(oo) described in 
text (dashed lines). (Bottom) Comparison with the marginal 
Fermi liquid theory - dashed lines. From [l| 

1 f u 

uom tr (uo) = - / dw'2u'm*(u'). (18) 
^ Jo 

The physical meaning of Eq. ([T6|) is the following: in op- 
tical measurements one photon with the energy uo is ab- 
sorbed and two excited particles (electron and hole) are 
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FIG. 5: Dependence of jl r (oo,T) on Soo in Bi2Sr2CaCu20s 
on different temperatures for €oo — 4 (solid lines) and Soo — 1 
(dashed lines). From [561 ]. 

created above and below the Fermi surface. If the elec- 
tron has energy uo' and the hole uo — a/, then they relax as 
quasi-particles with the renormalized Co. Since uo' takes 
values < uo' < uo then the optical relaxation uo tr {oo) is 
the energy- aver aged uo{uo) according to Eq. ([T6|) . The fac- 
tor 2 is due to the two quasi-particles, electron+hole. At 
finite T, the generalization reads [I?}, 03 

uotr{oo) = — / dw f \\ — tif{oo') — itlf{oo — uo')}2Co(uo'). 
u Jo 

(19) 

2. Inversion of the optical data and a 2 r (uo)F(uo) 
In principle, the transport spectral function 
a 2 r (uo)F(uo) can be precisely extracted from cr(cj), 
i.e. 7 tr (a;), only at T = K, which follows from Eq.( 

m 

1 d 2 

- • (20) 

However, real measurements are performed at finite T 
(and also at T > T c ) and the inversion procedure is in 
principle an ill-posed problem since a 2 r (uo)F(uo) is the 
deconvolution of the inhomogeneous Fredholm integral 
equation of the first kind with the temperature depen- 
dent Kernel K2(uo 1 v,T) in Eq. ([T2|) . An ill-posed mathe- 
matical problem, like this one, is very sensitive to input 
since experimental data contain less information than one 
needs. This can cause the fine structure of a 2 r (uo)F(uo) 
gets blurred in the extraction procedures and it can be 
temperature dependent even when the true a 2 r (uo)F(uo) 
is T independent. In the context of HTSC cuprates, this 
problem was first studied in [l8[, pjl with the follow- 
ing results: (1) the extracted shape of a 2 r (uo)F(uo) in 
Y Ba2CuzOj- x is not unique and it is temperature de- 
pendent, i.e. at higher T > T c the peak structure is 
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FIG. 6: Experimental (solid lines) and calculated (dashed 
lines) of R(uj) in optimally doped YBCO [57] at T=100, 200, 
300 K (from top to bottom). Inset: the two reconstructed 
at r (u)F(u) at T=100 K. The phonon density of states F{uS) 
- dotted line. From 



smeared and only a single peak (slightly shifted to higher 
uj) is present. For instance, the experimental data of 
R{uj) in YBCO were reproduced by two different spectral 
functions a 2 r (uj)F(uj), one with single peak and the other 
with three peaks structure as it is shown in Fig. [6] The 
similar situation is realized in optimally doped BISCO 
as it is seen in Fig. It is important to stress that the 
width of the extracted a 2 r (uj)F(uj) in both compounds 
coincide with the width of the phonon density of states 
F P h(v) EH, EH, 0; (2) the upper energy bound for 
a 2 r (uj)F(uj) can be extracted with certainty and it coin- 
cides approximately with the maximal phonon frequency 
in cuprates uj™* < 80 meV as it is seen in Figs. [6]|7l 

These results undoubtedly demonstrate the impor- 
tance of EPI in cuprates PH, pjl, [l|. We point out 
that the width of a 2 r (uj)F(uj) which is extracted from 
the optical measurements jl8| , pjl , [l| coincides with the 
width of the quasi-particle spectral function cx 2 {uj)F(uj) 
obtained in tunneling and ARPES spectra (which we 
shall discuss below), i.e. both functions are spread over 
the energy interval < uj < ^^)f x (^ 80 meV). Since 
in cuprates this interval coincides with the width in the 
phononic density of states F(uj) and since the maxima 
of a 2 (uj)F(uj) and F(uj) almost coincide, this is further 
strong evidence for the dominance of EPI. 

To this end, we would like to comment on two im- 
portant points. First, in some reports [59[,[l2[, [60[ it 
was assumed that a 2 r (uj)F(uj) of cuprates can be ex- 
tracted also in the superconducting state by using Eq. 
([20]) . However, Eq. ([20]) holds exclusively in the nor- 
mal state (at T=0) since a(uj) can be described by the 



FIG. 7: Experimental (solid line) and calculated (dashed line) 
of R(uj) in optimally doped BISCO [58] at T=100 K. Inset: 
the reconstructed o% r (oj)F(oj) - solid line. The phonon density 
of states F(uS) - dotted line. From [l8[ 



generalized (extended) Drude formula in Eq. JTTJ) only 
in the normal state. Such an approach apparently does 
not hold in the superconducting state since the dynami- 
cal conductivity depends not only on the electron-boson 
scattering but also on coherence factors and on the mo- 
mentum and energy dependent order parameter A(k, uj). 
In such a case it is unjustified to extract a 2 r (uj)F(uj) from 
Eq. ([20]) . Second, if R(uj) (and <j(uj)) in cuprates are due 
to some other bosonic scattering which is pronounced up 
to much higher energies uj c ^> Cc^ ax , this should be seen 
in the extracted spectral function a 2 r (uj)F(uj). Such an 
assumption is made, for instance, in the phenomenologi- 
cal spin-fluctuation approach p where it is assumed that 
a 2 (uj)F(uj) = g 2 flmx(oo) where lmx(uj) is extended up to 
the large energy cutoff uj c « 400 meV. This assumption 
is apparently in conflict with the above theoretical and 
experimental analysis which shows that solely EPI can 
describe R(uj) very well and that the contribution from 
higher energies uj ^> c<;£}f x must be small and therefore ir- 
relevant for pairing [l8| , [l9[ , [56[ . This is also confirmed 
by tunnelling measurements - see below. 

Despite the experimentally established facts for the im- 
portance of EPI, some reports appeared recently claiming 
that SFI dominates and that a 2 r (uj)F(uj) « g 2 ^lmx(^) 
where Imx(w) = / d 2 kx(k, uj) P3, 0- This claim is 
based on reanalyzing some old IR measurements [l2j , 
[§3]. The transport spectral function a 2 r (uj)F(uj) is ex- 
tracted in [l2[ by using the maximum entropy method in 
solving the Fredholm equation. However, in order to ex- 
clude negative values in the extracted a 2 r (uj)F(uj) they 
imply a biased condition that a 2 r {uj)F(uj) has a rather 
large tail at large energies - up to 400 meV. Then it is 
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not surprising at all that their extracted a 2 r {uo)F{uo) at 
large uo resembles qualitatively Imx(^) obtained by the 
magnetic neutron scattering on La2- x Sr x CuO^ [6l[. In 
other words one obtains as output what is assumed in 
the input. It turns out that even such a biased assump- 
tion in [l2[ by extracting a 2 r (uo)F(uo) does not reproduce 
the experimental curve Imx(w) [6l| in some important 
respects. (1 ) The relative heights of the two peaks in the 
extracted spectral function a 2 r (uo)F(uo) at lower temper- 
atures are opposite to that in Imx(w) [6l| - see Fig. 1 
in [l2[. (2) The strong temperature dependence of the 
extracted a 2 r (uo)F(uo), found in [12[, [60[ is in fact not 
the intrinsic property of the spectral function but it is 
due to the high sensitivity of the extraction procedure 
on temperature. As we already explained before, this is 
due to the ill-posed problem of solving the Fredholm in- 
tegral equation of the first kind with strong T-dependent 
kernel. (3) The extracted spectral weight a 2 r (uo)F(uo) 
in [l2[ has much smaller values at larger frequencies 
(uo > 100 meV) than it is the case for the measured 
ImxM, i.e. (I(uo > 100 meV)/I(u max )) <Clmx(^ > 100 
meV) /lmx(uo max ) - see Fig. 1 in [12j. In spite of the fact 
that the main weight of the extracted a 2 r (uo)F(uo) [l2[ 
lies in the range of phononic frequencies, < uo < uof^ ax 
it is not in agreement with that obtained in tunnelling 
and ARPES measurements. (4) To this end it is suspi- 
cious that the transport coupling constant X tr extracted 
in [l2[ is so large, i.e. X tr > 3 contrary to the previous 
findings that X tr < 1.5 [Hj], (H3|. Since in HTSC 
one has A > X tr this would probably give A « 6 that 
is not confirmed by other experiments. It is necessary 
to stress that the estimated X tr depends strongly on the 
value of plasma frequency, i.e. on uo 2 t - see Fig. [12] be- 
low, and it might be that for the latter the larger value 
is assumed in [12|. (5) The interpretation of a 2 r (uo)F(uo) 
in LSCO and BISCO solely in terms of Imx(^) is in con- 
tradiction with the magnetic neutron scattering in the 
optimally doped and slightly under doped YBCO [16j - 
that was discussed above, where in the former Imx(Q, uo) 
is small in the normal state - it is even below the exper- 
imental noise. This means that if the assumption that 
a 2 r (uo)F{uo) « g^lmx(uo) were correct then the contri- 
bution to lmx(uo) from the momenta < k < Q would 
be dominant and very detrimental for d-wave supercon- 
ductivity and T c would be rather low. Since the results 
of magnetic neutron scattering in YBCO [16[ are very 
convincing and trustful, then the conclusion is that the 
SFI coupling constant A s *^(~ g 2 ^) must be small, i.e. 

g sf < 0.2 eV and A s/ < 0.2. The latter is in accor- 
dance with other independent estimates (discussed also 
above) of A s/ (<C 2). 

Finally, we point out that very similar (to cuprates) 
properties, of <j(uo), R(uo) (and p(T) and electronic Ra- 
man spectra) were observed in 3D isotropic metallic ox- 
ides Lao.^Sro.sCoOs and Cclq^Stq^RuOs which are non- 
superconducting [62[ and in Ba\- x K x BiO^ which super- 
conducts at T c ~ 30 K at x = 0.4. This means that in 
all of them, the scattering mechanism might be of simi- 



lar origin. Since in these compounds there are no signs 
of antiferromagnetic fluctuations (which are present in 
cuprates), then EPI plays important role. 
3. Restricted optical sum-rule 

The restricted optical sum-rule was studied intensively 
in HTSC cuprates. It shows peculiarities not present in 
low-temperature superconductors. It turns out that the 
restricted spectral weight W(Vt Cl T) is strongly tempera- 
ture dependent in the normal and superconducting state, 
that was interpreted either to be due to EPI [20[, [2l| or 
to some non-phononic mechanisms [63|. In the follow- 
ing we demonstrate that the temperature dependence of 
W(Q C , T) = W(0) - pT 2 in the normal state can be ex- 
plained in a natural way by the T-dependence of the EPI 
transport relaxation rate 7^(0;, T) |20| , [2lj . Since the 
problem of the restricted sum-rule attracted much inter- 
est, it will be considered here in some details. In fact 
there are two kinds of sum rules related to cr(uo). The 
first one is the total sum rule which in the normal state 
reads 



Jo 



af {uo)duo = 



True" 
2m 



(21) 



while in the superconducting state it is given by the 
Tinkham-Ferrell- Glover (TFG) sum-rule 



a x (uo)duo = 
oX L 



af(oo)duo ■ 



(22) 



+0 



Here, n - the total electron density, e - the electron 
charge, m - the bare electron mass and Xl - the London 
penetration depth. The first (singular) term c 2 /8A| is 
due to the superconducting condensate which contributes 
o~f cond {oo) = (c 2 / 4X^)5(00). The total sum rule repre- 
sents the fundamental property of matter - the conserva- 
tion of the electron number, and to calculate it one should 
use the total Hamiltonian H tot = T e +Hi nt where all elec- 
trons, bands and their interactions Hint (Coulomb, EPI, 
with impurities, etc.) are accounted for. Here, T e is the 
kinetic energy of bare electrons 



/ a t p 2 
d 3 x2p (J (x)—ilj (J (x) 



2m 



y- 



4 



2m P CpffCp<T 



• (23) 



The partial sum rule is related to the energetics in 
the conduction (valence) band which is described by the 
Hamiltonian of the valence (band) electrons 



H 1} 



^ ^ GpC^p^-C^pcr -|- Vy^Coul- 



(24) 



p,cr 



It contains the band-energy with the dispersion e p and 
the effective Coulomb interaction of the valence electrons 
K,cw- In this case the partial sum-rule in the normal 
state reads [64] (for general form of e p ) 



^0 



a^ v (uo)duo = ^yYl 



(25) 
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where the number operator n V)P = ^ a c pa c pa ; l/m p = 
d 2 e p /dp 2 is the reciprocal mass and V is volume . In 
practice measurements are performed up to finite fre- 
quency and the integration over uo goes up to some cutoff 
frequency Q c (of the order of the band plasma frequency). 
In this case the restricted sum-rule has the form 



measure of the averaged kinetic energy. In the supercon- 
ducting state the partial sum-rule reads 



sxi 



+o 



af :V (uj)duj 



7re 2 a 2 



W(Sl c ,T)= / *?Ju)<Lj 



= -[K d + U(0)] 



ImU(u) 



duo. 



(26) 



where K d is the diamagnetic Kernel and II (lj) is the para- 
magnetic (current-current) response function - see more 
in [20[ , [2l| . In the case when the interband gap E g is the 
largest scale in the problem, i.e. when W b < ^ c < E g , in 
this region one has approximately ImII(a;) ~ and the 
limit Q c — > oo in Eq. ([26|) is justified. In that case one has 
11(0) « £ c (lmU(cj)/u)du which gives the approximate 
formula for W(Q C ,T) 



w(n c ,T) 



Jo 



(co)du; « -K d 



■ e 2 7r 



dp 2 



(27) 



where e p is the band-energy and n p = {n VyP ) is the quas- 
iparticle distribution function in the interacting system. 
Note that the right hand side of Eq. (j27|) does not depend 
on the cutoff energy Q c . So one should be careful not 
to interpret blindly the experimental result in cuprates 
by this formula and for that reason the best way is to 
calculate W(Vt Cl T) by using the exact result in Eq.([26|) 
which apparently depends on Q c . However, Eq. ([27|) is 
useful for appropriately chosen since it allows us to 
get semi-quantitative and qualitative results. In most 
papers related to the restricted sum-rule in HTSC, it 
was assumed, due to simplicity, the tight-binding model 
with nearest neighbors (n.n.) with the energy e p = 
— 2t(cosp x a + cosp y a) and l/m p = — 2ta 2 cosp x a. It is 
straightforward to show that in this case one has 



W(Sl c ,T) 



Jo ' 



v (u))dw 



ire 2 a 2 
2V 



(-T v ), 



(28) 



where (T v )h v — ^2 p ^p{ti v )h v is the averaged kinetic en- 
ergy of the band electrons and u p i jV is the (band) plasma 
frequency. In this approximation W(VL C1 T) is a direct 



(~T v )s. 



(29) 



In order to introduce the reader to the complexity of 
the problem of T-dependence of W(fi c ,T), let us con- 
sider the electronic system in the normal state and in 
absence of quasi-particle interaction. In that case one 
has n p = f p (/ p is the Fermi distribution function) 
and W n (fl c ,T) increases with the decrease temperature, 
i.e. W n (n c ,T) = W n (0) - P b T 2 where h ~ 1/W b . To 
this end, let us mention in advance that the experimen- 
tal value /3 exp is much larger than /3 b , i.e. /3 exp ^> f3 b 
thus telling us that the simple Sommerfeld-like smearing 
of f p by the temperature effects cannot explain the In- 
dependence of W(Vt c1 T) as it was put forward in some 
papers. We stress that the smearing of f p by temperature 
lowers the spectral weight compared to that at T = 
i.e. W n (Q c ,T) < W n (Q c ,0). In that respect it is not 
surprising at all that there is a lowering of W S (Q C ,T) 
in the BCS superconducting state, W^ cs (ft c ,T = 0) < 
W n (fl c ,0) since f p is smeared due to the appearance of 
the superconducting gap, 2f p = 1 — (£ p /E p )th(E p /2T), 

E p = y^£ P + A 2 , ^ p = e p — /i, and the maximal decrease 

of W s (n c ,T) is atT = 0. 

Let us enumerate and analyze the main experimental 
results in cuprates. 1. In the normal state (T > T c ) of 
most cuprates, one has W n (Q c ,T) = W n (0) — j3 ex T 2 with 
/^exp ^ ftbi i- e - W n (Q c ,T) is increasing by decreasing 
T, even at T below the opening of the pseudogap. The 
change of W n (^ c , T) from room temperature down to T c 
is no more than 5 %. 2. In the superconducting state (T < 
T c ) of some underdoped and optimally doped Bi-2212 
compounds [66[, [65[ (and underdoped Bi-2212 films [72| ) 
there is an effective increase of W S (Q C ,T) with respect 
to that in the normal state, i.e. W 3 (Cl c ,T) > W n (Vt c1 T) 
for T < T c . This is non-BCS behavior shown in Fig. [U 

In some optimally doped and in most overdoped 
cuprates, there is decreasing of W 3 (£l c ,T) at T < T c 
which is the BCS-like behavior [67] as it is seen in Fig. [9] 

We stress that the non-BCS behavior of W S (Q C ,T) 
for underdoped and optimally doped systems was ob- 
tained by assuming that ^ c ~ (1 — 1.2) eV. How- 
ever, in Ref. [55[ these results have been questioned 
and the conventional BCS-like behavior was observed 
(W s (Sl c ,T) < W n (Sl c ,T)) in the optimally doped YBCO 
and slightly underdoped Bi-2212 by using larger cutoff 
energy Q c = 1.5 eV. Although the results obtained in 
[55l | looks very trustfully, it is fair to say that the issue of 
the reduced spectral weight in the superconducting state 
of cuprates is still unsettled and under dispute. In over- 
doped Bi-2212 films, the BCS-like behavior W s (Vt c ,T) < 
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FIG. 8: Measured spectral weight W s {Vt c , T) for uj c w 1.25eV 
in two underdoped Bi2212 (with T c = 88 K and T c = 66 K). 
From [66j ]. 
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W n (fi c , T) was observed, while in LS CO it was found that 
W s (n c ,T) * const, i.e. W S (^ C ,T < T c ) w W n (ft c ,T c ). 

How to explain the strong temperature dependence of 
W{Vt c ,T) in the normal and superconducting state? In 
Hfli, HH] it was shown that the EPI relaxation j ep (T) 
plays the main role in the T-dependence of W(Q C ,T). 
The main theoretical results of [20], [2l| are the follow- 
ing. (1) The calculations based on the exact formula in 
Eq. ([27j) give that for Q c ^> Qd, the difference in spectral 
weights of the normal and superconducting state is small, 
i.e. W n (n c ,T) « W s (n c ,T) (in the following we call it 
W) since W n (Q c ,T) - W s (n c ,T) ~ A 2 /n 2 c . In the case 
of large Q c based on the approximate formula Eq. ([27|h 
one obtains 

In the case of EPI, one has 7 = j ep (T) + 7* m ^ where 
7 ep ( T ) = dza 2 (z)F(z)coth(z/2T). It turns out 
that for a 2 (u)F(u) shown in Fig. [3l one obtains: 
(i) >y ep ( T ) ~ ^ 2 in the temperature interval 100 



FIG. 9: (Top) Spectral weight W n (ftc,T) of the overdoped 
Bi2212 for Q c = leV. Closed symbols - normal state. Open 
symbols - superconducting state. (Bottom) Change of the 
kinetic energy AE k i n = E k i n ,s - E k in,N in meV per Cu site 
vs the charge p per Cu with respect to the optimal value p op t- 
From |3. 



K < T < 200 K as it is seEn in Fig. [TUl HJ, jp; 
(ii) the second term in Eq. ([30|) is much larger than 
the last one (the Sommerfeld-like term). For the EPI 
coupling constant = 1.5 one obtains rather good 
agreement with experiments. At lower temperatures, 
7 ep (T) deviates from the T 2 behavior and the deviation 
depends on the structure of the spectrum in a 2 {uo)F(uo). 
It is seen in Fig. [10] that for a softer Einstein spectrum 
(with tt E = 200 K), W(Q C ,T) lies above the curve 
with the T 2 asymptotic, while the one with a harder 
phononic spectrum (with Qe = 400 K) lies below it. 
This result means that different behavior of W(fi c ,T) 
in the superconducting state of cuprates for different 
doping might be simply related to different contributions 
of low and high frequency phonons. We stress that such 
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FIG. 10: Spectral weight W(Q C , T) for Einstein phonons with 
Q E = 200 K (full triangles) and Qe = 400 K (open circles, 
left axis). Dashed lines is T 2 asymptotic. From [211 ]. 



a behavior of W(f2 c ,T) was observed in experiments 
[66[, [65[, [55[ and the above analysis tells us that the 
theory based on EPI explains in a consistent way the 
strange temperature behavior of W(Q C ,T) above and 
below T c and that there is no need to invoke exotic 
scattering mechanisms. 

4. Resistivity p(T) 

The temperature dependence of the in-plane resistivity 
p ab (T) in cuprates is a direct consequence of the quasi-2D 
motion of quasi-particles and of the inelastic scattering 
which they experience. At present, there is no consensus 
on the origin of the linear temperature dependence of the 
in-plane resistivity p ab (T) in the normal state and there is 
rather widespread believe that it can not be due to EPI. 
The inadequacy of this belief was already demonstrated 
by analyzing the dynamic conductivity <j(oo) which is suc- 
cessfully explained by EPI. Since p(T) = l/a(u = 0) 



p(T) = % ltr {T) 



cut 



J imp 



(31) 



7trC0 



duj- 



a 2 tr (u)F(u). (32) 



It is quite natural that in some temperature region, p(T) 
in cuprates can be explained by EPI as it is shown in Fig. 
QU It turns out that % r (T) ~ T for T > aG D , a < 1 
depending on the shape of a^ r {uo)F{uo). In case of the 
Debye spectrum, it is realized for T > Od/5 i-e. 



P(T) 



87r K h4 =pT - 



(33) 
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FIG. 11: (a) Calculated resistivity p(T) for the EPI spectral 
function al r {uj)F{uj) in [68]. (b) Measured resistivity in a(x)- 
and b(y)-crystal direction of YBCO [70|] and calculated Bloch- 
Griineisen curve for X ep = 1, [711 ]. 



There is an experimental constraint on \ tr , i.e. 
X tr w 0.2bu;l l (eV)p , (fjincm/K), 



(34) 



which imposes a limit on it. For instance, for co p i « (2—3) 
eV [62[ and p' ~ 0.6 in the oriented YBCO films and 
p r w 0.3 — 0.4 in single crystals of BSCO, one obtains 
X tr ~ 0.4 — 1.2. In case of YBCO single crystals, there 
is a pronounced anisotropy in p ab (T) [70[ which gives 
p' x (T) = 0.6pQcm/K and p' y (T) = 0.25/i^cm/K. Ac- 
cording to Eq. ([34|) , one obtains Xtr pi) which is shown 
in FigH2j where the plasma frequency uo p i which enters 
Eqs. (|31ti33|) can be calculated by LDA and also extracted 
from the width (~ lu*^ of the Drude peak at small fre- 
quencies, where uo p i = y/e^u*^ 

We shall argue below that from tunnelling experiments 
[23[ -[2^ one obtains in the framework of the Eliashberg 
theory that the EPI coupling constant is large X ep ^2 — 3 
which implies that X tr ~ (X/S), i.e. EPI is reduced in 
transport properties due to some reasons that shall be 
discussed in Part II. Such a large reduction of X tr cannot 
be obtained within the LDA band structure calculations 
which means that X ep and Xtr contain renormalization 
which do not enter in the LDA theory. In Part II we 
shall argue that the strong suppression of X tr may have 
its origin in strong electronic correlations and the long- 
range Madelung energy [45[, [73j]. 

4. Femtosecond time- resolved optical spectroscopy 
The femtosecond time-resolved optical spectroscopy 
(FTROS) has been developed intensively in the last cou- 
ple of years and applied successfully to HTSC cuprates. 
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FIG. 12: Transport EPI spectral function coupling constant 
in YBCO as a function of plasma frequency u p as derived from 
the experimental slope of resistivity p'(T) in Eq. (|34|) . A^ for 
p' x (T) = 0.6pQcm/K and \ y for p' y {T) = 0.25/xQcm/if 0. 
Squares are LDA values [691 ]. 



In this method a femtosecond (lfs = 10 -15 sec) laser 
pump excites in materials electron-hole pairs via inter- 
band transitions. These hot carriers release their en- 
ergy via electron-electron (with the relaxation time r ee ) 
and electron-phonon scattering reaching states near the 
Fermi energy within 10 — 100 fs - see [74]. The typical 
energy density of the laser pump pulses with the wave- 
length A « 810 nm (1.5 eV) was around F ~ lpj/cm 2 
(the excitation fluenc F) which produces approximately 
3 x 10 10 carriers per puis (by assuming that each photon 
produces huu/A carriers, A is the superconducting gap). 
By measuring photoinduced changes of the reflectivity in 
time, i.e. AR(t)/Ro, one can extract information on the 
further relaxation dynamics of the low-laying electronic 
excitations. Since AR(t) relax to equilibrium the fit with 
exponential functions is used 



AR(t) 
Ro 



fit) 



Ae t a + Be~ 



(35) 



where f(t) = H(t)[l - exp{-£/r ee }] (H(t) is the Heavy- 
side function) describes the finite rise-time. The parame- 
ters A, B depends on the fluenc F. This method was used 
in studying the superconucting phase of La2- x Sr x CuO±, 
with x = 0.1, 0.15 and T c = 30 K and 38 K respec- 
tively UK. In that case the signal A ^ for T<T C and 
A = for T > T c , while the signal B was present also 
at T > T c . It turns out that the signal A is related 
to the quasi-particle recombination across the supercon- 
ducting gap A(T) and has a relaxation time of the order 
r A > 10 ps at T = 4.5 K. At the so called thresh- 
old fluenc {F T = 4.2 ± 1.7 pJ/cm 2 for x = 0.1 and 
Ft = 5.8 ± 2.3 /a J /cm 2 for x = 0.15) occurs the vapor- 
ization (destroying) of the superconducting phase, where 
the parameter A saturates. This vaporization process 



takes place at the time scala r r « 0.8 ps. The exter- 
nal fluenc is distributed in the sample over the excitation 
volume which is proportional to the optical penetration 
depth \ op (~ 150 nm at A ~ 810 nm) of the pump. The 
energy densities stored in the excitation volume at the 
vaporization threshold for x = 0.1 and x = 0.15 are 
U p = F T /X p = 2.0 ± 0.8 K/Cu and 2.6 ± 1.0 K/Cu, re- 
spectively. The important fact is that U p is much larger 
than the superconducting condensation energy which is 
U cond « 0.12 K/Cu for x = 0.1 and U cond « 0.3 K/Cu 
for x = 0.15, i.e. U p ^> U con d- This means that 
the energy difference U p — U con d must be stored else- 
where on the time scale r r . The only present reservoir 
which can absorb the difference in energy are the bosonic 
baths of phonons and spin fluctuations. The energy re- 
quired to heat the spin reservoir from T = 4.5 K to T c is 
U 8f = C 8f (T)dT. The measured C sf (T) in La 2 Cu0 4 
|22j gives very small value U s f « 0.01 K. In the case of 

the phonon reservoir on obtains U p h = J^ c C p h(T)dT = 9 
K/Cu for x = 0.1 and 28 K/Cu for x = 0.15. Since 
U s f <C Up — U con d the spin reservoir cannot absorb the 
rest energy U p — U con d' The situation is opposite with 
phonons since U p h ^>U p — U con d and phonon can absorb 
the rest energy in the excitation volume. The complete 
vaporization dynamics can be described in the framework 
of the Rothwarf- Taylor model which describes approach- 
ing of electrons and phonons to quasi-equilibrium on the 
time scale of 1 ps [75]. We shall not go into details but 
only summarize, that only phonon- mediated vaporiza- 
tion is consistent with the experiments, thus ruling out 
spin-mediated quasi-particle recombination and pairing 
in HTSC cuprates. This is additional proof for the inef- 
fectivness of the SFI scattering in cuprates. 

In conclusion, optics and resistivity measurements in 
normal state of cuprates are much more in favor of EPI 
than against it. However, some intriguing questions still 
remain to be answered: (i) what are the values of \ tr 
and uo p i] (ii) what is the reason that X tr <C A is realized 
in cuprates; (iii) what is the role of Coulomb scattering 
in a(u) and p(T). Later on we shall argue that ARPES 
measurements in cuprates give evidence for a contribu- 
tion of Coulomb scattering at higher frequencies, where 
7(0;) ~ 7 + A c cj for uo > c^^ ax with A c w 0.4. So, despite 
the fact that EPI is suppressed in transport properties 
it can be sufficiently strong in the self-energy in some 
frequency and temperature range. 



C. ARPES and the EPI self-energy 

ARPES is nowadays a leading spectroscopy method in 
the solid state physics [ll|. It provides direct informa- 
tion to the one-electron removal spectrum in a complex 
many system. The method involves shining light (pho- 
tons) with energies between 5 — 1000 eV on the sam- 
ple and by detecting momentum (k) - and energy(o;)- 
distribution of the outgoing electrons. The resolution 
of ARPES has been significantly increased in the last 
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decade with the energy resolution of AE ^1 — 2 meV 
(for photon energies ~ 20 eV) and angular resolution 
of AO < 0.2°. The ARPES method is surface sensitive 
technique, since the average escape depth (l esc ) of the 
outgoing electrons is of the order of l esc ~ 10 A, depend- 
ing on the energy of incoming photons. Therefore, very 
good surfaces are needed in order that the results be rep- 
resentative for bulk samples. The most reliable studies 
were done on the bilayer Bi 2 Sr2CaCu20s (Bi2212) and 
its single layer counterpart B%2Sr2CuO% (£>i2201), since 
these materials contain weakly coupled BiO planes with 
the longest inter-plane separation in the cuprates. This 
results in a natural cleavage plane making these materials 
superior to others in ARPES experiments. After a drastic 
improvement of sample quality in other families of HTSC 
materials, the ARPES technique has became a central 
method in theoretical considerations. Potentially, it gives 
valuable information on the quasi-particle Green's func- 
tion, i.e. on the quasi-particle spectrum and life-time ef- 
fects. The ARPES can indirectly give information on the 
momentum and energy dependence of the pairing poten- 
tial. Furthermore, the electronic spectrum of the (above 
mentioned) cuprates is highly quasi-2D which allows an 
unambiguous determination of the initial state momen- 
tum from the measured final state momentum, since the 
component parallel to the surface is conserved in photoe- 
mission. In this case, the ARPES probes (under some 
favorable conditions) directly the single particle spectral 
function A(k, uj). In the following we discuss only those 
ARPES experiments which give evidence for the impor- 
tance of the EPI in cuprates - see more in [111 ]. 

The photoemission measures a nonlinear response 
function of the electron system, and under some con- 
ditions it is analyzed in the so-called three-step model, 
where the total photoemission intensity It t(k, uj) w 
/ • I2 • I3 is the product of three independent terms: 
(1) I - describes optical excitation of the electron in 
the bulk; (2) I2 - describes the scattering probabil- 
ity of the travelling electrons; (2) I3 - the transmis- 
sion probability through the surface potential barrier. 
The central quantity in the three-step model is /(k, uj) 
and it turns out that it can be written in the form 
(for k = k,|) [11] J(k,cj) ~ I (Kv)f(uj)A(k,uj) with 
Iq (k, v) ~| (ipf I pA |^) 1 2 and the quasi-particle spec- 
tral function A(k, uj) = — ImG(k, uj)/tt 

1 IraX(k,o;) 

[ ,<J) ~ ~n[uj-£(k) -ReX(k,o;)] 2 + JraX 2 (k,o;)' 

(36) 

Here, {ip ^ | pA |^) is the dipole matrix element which 
depends on k, polarization and energy v of the incoming 
photons. The knowledge of the matrix element is of a 
great importance and its calculation from first principles 
was done carefully in [76]. f(ui) is the Fermi function, 
G and E = Re E + i Im E are the quasi-particle Green's 
function and the self-energy, respectively. 

We summarize and comment here on some important 
ARPES results which were obtained recently and which 



confirm the existence of the Fermi surface and impor- 
tance of EPI in quasi-particle scattering [ll| . 
ARPES in the normal state 

(Nl) There is a well defined Fermi surface in the metal- 
lic state with the topology predicted by the LDA. How- 
ever, the bands are narrower than LDA predicts which 
points to a strong quasi-particle renormalization. (N2) 
The spectral lines are broad with |ImX(k, uj) |~ uj (or 
~ T for T > uj) which tells us that the quasi-particle 
liquid is a non-canonical Fermi liquid. (N3) There is a 
bilayer band splitting in Bi22\2 (at least in the over- 
doped state). The previous experiments did not show 
this splitting and served for various speculations on some 
exotic non- Fermi liquid scenarios. (N4) At temperatures 
T c < T < T* and in the under-doped cuprates there is a 
d-wave like pseudogap A p ^(k) ~ A p ^o(cos k x — cos k y ) in 
the quasi-particle spectrum where A PQi o increases by low- 
ering doping. We stress that the pseudogap phenomenon 
is not well understood at present and we shall discuss this 
problem in Part II. Its origin can be due to a precursor 
superconductivity or due to a competing order, such as 
spin- or charge-density wave or something similar. (N5) 
The ARPES self-energy gives clear evidence that EPI 
interaction is rather strong. For instance, at T > T c 
there are kinks in the quasi-particle dispersion o;(£ k ) in 
the nodal direction (along the (0, 0) — (tt, tt) line) at the 
characteristic phonon energy uj^ ~ (60 — 70) meV [48^ 
see Fig. [13j and near the anti-nodal point (tt, 0) at 40 
me V [77] - see Fig. [131 

That these kinks exist also above T c excludes the sce- 
nario with the magnetic resonance peak in Jmx s (Q,w). 
Since the magnetic neutron scattering give small SFI cou- 
pling constant X s ^ < 0.3. The kinks cannot be due to 
SFI as we already discussed above. (iV6) The position of 
the nodal kink is practically doping independent which 
points towards phonons as the scattering (gluing) boson. 
(N7) The quasi-particles (holes) couple practically to the 
whole spectrum of phonons since at least three group of 
phonons were extracted from the ARPES effective self- 
energy in La 2 - x Sr x Cu04 - Fig. [H 

This result is in a qualitative agreement with numer- 
ous tunnelling measurements [23[-[28| which apparently 
demonstrate that the broad spectrum of phonons cou- 
ples with holes without preferring any particular phonons 
- see discussion below. (NS) Recent ARPES measure- 
ments in B2212 [49| show very different slope duj/d^ of 
the quasi-particle energy u;(£ k ) at very small | £ k |<C oj p h 
and large energies | £ k |^> uj p h - see Fig. [151 The theoret- 
ical analysis [79[ of these results gives the total coupling 
constant A > 3, the EPI one X ep > 2 while the Coulomb 
scattering (SFI is a part of it) is A c « 1 [79] - see Fig. [15l 
To this end let us mention some confusion related to 
the value of the EPI coupling constant extracted from 
ARPES. Namely, in [n|, Hflf, H3 the EPI self-energy 
was obtained by subtracting the high energy slope of the 
quasi-particle spectrum uj(^ k ) at uj ~ 0.3 eV. The latter 
is apparently due to the Coulomb interaction. Although 
the position of the low-energy kink is not affected by this 
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FIG. 14: (a) Effective real self-energy for the non- 
superconducting La2-xSr x Cu04, x = 0.03. Extracted 
a 2 eff {uj)F(uj) is in the inset, (b) Top: the total MDC width - 
open circles. Bottom: the EPI contribution shows saturation, 
impurity contribution - dotted black line. The residual part 



is growing 



From [78[. 



stant. The quasi-particle spectrum uo{k) is determined 
from the condition 



■f(k) -Re[E e ^(k,cj) -f E c (k,c 



0, 



(37) 



where £(k) is the bare band structure energy. At low 
energies uo < cc^ ax <C uo Cl Eq. (|37j) can be rewritten in the 
form 



FIG. 13: (Top) Quasi-particle dispersion of BiTlYl, 5*2201 
and LSCO along the nodal direction, plotted vs the momen- 
tum k for (a) — (c) different doping, and (d) — (e) different T; 
black arrows indicate the kink energy; the red arrow indicates 
the energy of the q = (tt, 0) oxygen stretching phonon mode; 
inset of (e)- T-dependent E ; for optimally doped Bi2212; (f) - 
doping dependence of the effective coupling constant A' along 
(0,0) - (7r,7r) for the different HTSC oxides. From Ref. 0. 
( Bottom) Quasi-particle dispersion E(k) in the normal state 
(al, bl, c), at 107 K and 115 K, along various directions cj) 
around the anti-nodal point. The kink at E — AQmeV is 
shown by the horizontal arrow. (a2 and b2) is E(k) in the su- 
perconducting state at 10 K with the shifted kink to 70meV. 
(d) kink positions as a function of <j) m the anti-nodal region. 
From Ref. Fil 



procedure (if c^ ax <C uo c ), the above (subtraction) proce- 
dure gives in fact an effective EPI self energy E^.(k, uo) 



and 



The latter 



the coupling constant A^ p e jj(k) only, 
is smaller than the real EPI coupling constant A^ p (k). 
The total self-energy is E(k,a;) = E ep (k,^) + E c (k,cj) 
where E c is the contribution due to the Coulomb inter- 
action. At very low energies uo <C uo c one has usually 
E c (k, uo) = — A^(k)u;, where uo c (~ 1 eV) is the character- 
istic Coulomb energies and X c z the Coulomb coupling con- 



a;-r n (k)-ReEf // (k,^) = 0, 



r n (k) = [i+A^(k)]- i e(k) 



ReZf ff (k,uo) 



Re^ ff (^uo) 
1 + A*(k) 



J ph 



(38) 



(39) 



(40) 



one has 



Since at very low energies uo <C 
ReE e ^(k,cj) = -Xl p (k)u and ReE^ / (k,^) = 
— A^ p e jj(k)cj, then the real coupling constant is related 
to the effective one by 

Af (k) = [1 + K(k)]\Z f ,(k) > \Zf f (k). 



At higher energies uo™ h ax < uo < uo c , which are less 
important for pairing, the EPI effects are suppressed 
and E ep (k, uo) stops growing, one has ReE(k, uo) ~ 
ReE e ^(k,cj) - X c z (k)uo. The measured ReE exp (k,cj) at 
T = 10 K near and slightly away from the nodal point 
in the optimally doped Bi2212 with T c = 91 K [Hi is 
shown in Fig. [15l 

It is seen that i?eE exp (k, uo) has £wo kinks - the first 



one at low energy uo\ 



high 
J ph 



50 — 70 meV which 



is most probably of the phononic origin [11|, [80[, [81|, 
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FIG. 15: Fig.4b from ReZ(u) measured in Bi2212 (thin 
line) and model ReY,(uj) (bold line) obtained in [83]. The 
three thin lines (Ai, A2, A3) are the slopes of ReY>(oj) in differ- 
ent energy regions - see the text. 
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FIG. 16: (Left) The ARPES spectrum of undoped 
CaiCuOiCli at k = (7r/2, 7r/2). Gaussian shape - solid line, 
Lorentzian shape - dashed line. (Right) Dispersion of the po- 
laronic band - A and of the quasi-particle band - B along the 
nodal direction. Horizontal lines are the chemical potentials 
for a large number of samples. From [84[ 



while the second kink at higher energy 002 ~ uj c « 350 
meV is probably due to the Coulomb interaction. How- 
ever, the important results in Ref. [83] is that the slopes 
of ReS exp (k, uj) at low (uj < u^ h ) and high energies 

(ujp^ h < oj < oj c ) are different. The low-energy and 
high-energy slope near the nodal point are depicted and 
shown in Fig. fT5l schematically (thin lines). From Fig. [15] 
it is obvious that EPI prevails at low energies uj < . 
More precisely digitalization of Re £ exp (k, uj) in the in- 
terval < uj < 0.4 eV gives the Coulomb coupling 
X c z ~ 1.1 while the same procedure at 20 meV ~ uj 1 °™ < 
uj < uJph h ~ 50 — 70meV gives the total coupling con- 
stant (A2 =)X Z = X ep + X c z « 3.2 and the EPI coupling 
constant Xf (= Xf Mgh ) « 2.1 > 2A£ e// (k), i.e. the 
EPI coupling is at least twice larger than in the previ- 
ous analysis of ARPES results. This estimation tells us 
that at (and near) the nodal point, the EPI interaction 
dominates in the quasi-particle scattering at low ener- 
gies since X ep {^ 2.1) w 2A^ > 2X S J , while at large en- 
ergies (compared to uj p h), the Coulomb interaction with 
X c z w 1.1 dominates. We point out that EPI near the anti- 
nodal point can be even larger than in the nodal point, 
mostly due to the higher density of states near the anti- 
nodal point. (NS) Recent ARPES spectra in the opti- 
mally doped Bi2212 near the nodal and anti- nodal point 
(8lj show a pronounced isotope effect in Re S exp (k, uj) : 
thus pointing to the important role of EPI - see more 
in the part related to the isotope effect. The isotope 
effect in ReS(k, uj) can be well described in the frame- 
work of the Migdal-Eliashberg theory for EPI [82] as it 
will be discussed in Part II. (N9) ARPES experiments 
on Ca2Cu02Cl2 give strong evidence for the formation 



of small polarons in undoped cuprates which can be only 
due to phonons and strong EPI, while by doping quasi- 
particles app ear and there are no small polarons [Uj]. 
Thus in [84], it a broad peak (around —0.8 eV) is ob- 
served at the top of the band (k = (7r/2,7r/2)) with the 
dispersion similar to that predicted by the t — J model - 
see Fig. [TBI 

However, the peak in Fig. [16] is of Gaussian shape and 
can be described only by coupling to bosons, i.e. this 
peak is a boson side band - see more in [5] and refer- 
ences therein. The theory based on the t-J model (in the 
antiferromagnetic state of the undoped compound) by 
including coupling to several (half-breathing, apical oxy- 
gen, low-lying) phonons, which is given in [85], explains 
successfully this broad peak of the boson side band by 
the formation of small polarons due to the EPI coupling 
(X ep « 1.2). Note that this X ep is for the polaron at the 
bottom of the band while in the case whwn the Fermi sur- 
face exists this coupling is even larger [85[ . In [85[ it was 
stressed that even when the electron-magnon interaction 
is stronger than the EPI one, the polarons are formed due 
to EPI. The latter involves excitation of many phonons 
at the lattice site (where the hole is seating), while it is 
possible to excite only one magnon at the given site. 

(A/10) Recent soft x-ray ARPES measurements 
on the electron doped HTSC Ndi, 85 Ce^^CuO^ [861 ]. 
and Sm i2 - x )Ce x Cu0 4 (x = 0.1, 0.15, 0.18), 
CcoasCuOa, Eui.sbCeo^CuO^ |87] show kink at 
energies 50 — 70 meV in the quasi-particle dispersion re- 
lation along both nodal and antinodal directions as it is 
shown in Fig. [TT1 

It is seen from this figure that the effective EPI cou- 
pling constant X ep ^A< X ep ) is isotropic and X e e p ^ « 0.8— 1. 
The kink in the electron-doped cuprates is due solely to 
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FIG. 17: NCCO electron-doped: (a) JmE(w) measured in 
the nodal point. Curves are offsets by 50 meV for clarity. 
The change of the slope in the last bottom curve is at the 
phonon energy, (b) JmE(w) for the antinodal direction with 
30 meV offset, (c) Experimental phonon dispersion of the 
bond stretching modes, (d) Estimated X e e p ^ from 
From 0. 



EPI and in that respect the situation is similar to the 
one in the hole-doped cuparates. 

ARPES results in the superconducting state 
(SI) There is an anisotropic superconducting gap in 
most HTSC compounds [ll|, which is predominately 
d-wave like, i.e. A sc (k) ~ Aofcos/c^ — cosk y ) with 
2A /T C « 5 - 6. (52) The kink at (60 - 70) meV in 
the quasi-particle energy around the nodal point is not- 
shifted in the superconducting state while the antinodal 



kink at cc/^ ~ 40 meV is shifted in the superconducting 

state by A (= (25 - 30)meV), i.e. uoffi -> + A = 
(65 — 70)meV [11]. To remind the reader, in the standard 
Eliashberg theory the kink in the normal state at uo = uo p h 
should be shifted in the superconducting state to u) p h+ Ao 
at any point at the Fermi surface. This puzzling result 
might be a smoking gun result since it makes a constraint 
on the quasi-particle interaction in cuprates. Until now 
there is only one plausible explanation [88] of this shift- 
non-shift puzzle which is based on an assumption of the 
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FIG. 18: Crystal structure of Ba2Ca 3 Cu40 8 (08F 1 -s)2. 
There are four CuOi layers in a unit cell with the outer having 
apical F atoms. From |89l ]. 



forward scattering peak (FSP) in EPI - see more in Part 
II. The FSP in EPI means that electrons scatter into a 
narrow region around the starting point in the k-space, 
so that at the most part of the Fermi surface, there is 
weak (or no) mixing of states with different signs of the 
order parameter A(k). (S3) The recent ARPES spectra 
[89[ on an undoped single crystalline 4-layered cuprate 
with apical fluorine (F), Ba2Ca^Cu^O^F2 (F0234) gives 
strong evidence against SFI - see Fig. [TBI 

Namely, F0234 is not a Mott insulator - as expected 
from valence charge counting which puts Cu valence as 
2 + , but it is a superconductor with T c = 60 K. More- 
over, the ARPES data [89[ reveal at least two metallic 
Fermi-surface sheets with corresponding volumes equally 
below and above half-filling - see Fig. [191 

One of the Fermi-surfaces is due to the electron-like 
(TV) band (with 20 ± 6% electron-doping) and the other 
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FIG. 19: (a) Fermi surface (FS) contours from two samples 
of F0234. n - electron-like; h - hole-like. Bold arrow is (7r,7r) 
scattering vector. Angle 6 defines the horizontal axis in (b). 
(b) Leading gap edge along k-space angle from the two FS 
contours. From [891 ]. 



one due to the hole-like (P) band (with 20 ± 8 % hole- 
doping) and their split along the nodal direction is sig- 
nificant and cannot be explained by the LDA (or DFA) 
method [90[. This electron and hole self-doping of in- 
ner and outer layers is in an appreciable contrast to 
other multilayered cuprates where there is only hole self- 
doping. For instance, in HgBa2Ca n Cu n+ i02 n +2 (n = 
2,3) and (Cu,C)Ba2Ca n Cu n+ i03 n+ 2 (n = 2,3,4), the 
inner Cu02 layers are less hole-doped than outer layers. 
It turns out, unexpectedly, that the superconducting gap 
on the Af-band Fermi-surface is significantly larger than 
on the P-one, i.e. their ratio is anomalous (Ajv/Ap) « 2 
and A 7v is an order of magnitude larger than in the 
electron-doped cuprate A d2- x C -e x CuO '4. Furthermore, 
the Af-band Fermi-surface is rather far from the antin- 
odal point at (7r/2,0). This is an extremely important 
result which means that the antiferromagnetic spin fluc- 
tuations with the AF wave- vector Q = (7r/2,7r/2), as 
well as the van Hove singularity, are not important for 
the pairing in the A'-band. To remind the reader, the SFI 
scenario assumes that the pairing is due to spin fluctu- 
ations with the wave- vector Q which connects two anti- 
nodal points which are near the van Hove singularity at 
the hole-surface (at (7r/2,0) and (0,7r/2)) giving rise to 
large density of states. This is apparently not the case 
for the A'-band Fermi-surface - see Fig. [T9j The ARPES 
data give further that there is a kink at ~ 85 meV in the 
quasi-particle dispersion of both bands, while the kink in 
the A'-band is stronger than that in the P-band. This 
result, together with the anomalous ratio (Ajv/Ap) « 2, 
strongly disfavors SFI as a pairing mechanism. 

(S4) Despite the presence of significant elastic scat- 
tering in optimally doped Bi-2212, there are dramatic 
sharpening of the spectral function near the anti-nodal 
point (7r,0) at T < T c [9l| . This can be explained by 



assuming that the small q-scattering (the forward scat- 
tering peak) dominates in the elastic impurity scattering 
[92| | . [93^ . As a result, one finds that the impurity scatter- 
ing rate in the superconducting state is almost zero, i.e. 

7im P ( k ^) = 7n( k ^) +7a( k ^) = for | | < A for 
any kind of pairing (s- p- d-wave etc.) since the normal 
(7 n ) and the anomalous (7 a ) scattering rates compensate 
each other - the collapse of the elastic scattering rate. 
This result is a consequence of the Anderson-like theo- 
rem for unconventional superconductors which is due to 
the dominance of the small q-scattering [92j]-[93j]. In such 
a case d-wave pairing is weakly unaffected by impurities 
- there is small reduction in T c [93[ , [94| . The physics be- 
hind this result is rather simple. The small q-scattering 
(forward scattering) means that electrons scatter into a 
small region in the k-space, so that at the most part of 
the Fermi surface there is no mixing of states with dif- 
ferent signs of the order parameter A(k), and the detri- 
mental effect of impurities is reduced. For states near 
the nodal points, there is mixing but since A(k) is small 
in this region, there is only small reduction in T c [42I ]. 
This result points to the importance of strong correla- 
tions in the renormalization of the impurity scattering - 
see discussion in Part II. 

In conclusion, in order to explain ARPES results in 
cuprates it is necessary to take into account: (1) EPI 
interaction since it dominates in the quasiparticle scat- 
tering in the energy region responsible for pairing; (2) 
effects of elastic nonmagnetic impurities with FSP; (3) 
the Coulomb interaction which dominates at higher en- 
ergies uo > uo p h- In this respect, the presence of ARPES 
kinks and the knee-like shape of the spectral width are 
serious constraints for the pairing theory. 



D. Tunnelling spectroscopy and spectral function 

By measuring current-voltage I — V characteristics 
in NIS (normal metal-insulator-superconductor) tun- 
nelling junctions with large tunnelling barrier - see 
more below, one obtains from tunnelling conductance 
Gns(V) = dl/dV the so called tunnelling density of 
states in superconductors Nt(uj). Moreover, by mea- 
suring of Gats' (V) at voltages eV > A it is possible to 
determine the Eliashberg spectral function a 2 F(u) and 
finally to confirm (definitely) the phonon mechanism of 
pairing in LTSC materials, except heavy fermions. Four 
tunnelling techniques were used in the study of HTSC 
cuprates: (1) vacuum tunneling by using the STM tech- 
nique - scanning tunnelling microscope; (2) point-contact 
tunnelling; (3) break-junction tunnelling; (4) planar- 
junction tunnelling. Each of these techniques has some 
advantages although in principle the most potential one 
is the STM technique [95[. It should be stressed that 
there are still difficulties in understanding tunnelling ex- 
periments in HTSC cuprates because of non-ideal tun- 
nelling behavior of contacts [95j]. Since tunnelling mea- 
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surements probe a surface region of the order of supercon- 
ducting coherence length £ , then this kind of measure- 
ments in HTSC materials with small coherence length 
£o (£ab ^ 20 A in the a — b plane and £ c ~ 1 — 3 A along 
the c — axis) depends strongly on the surface quality and 
sample preparation. Nowadays, many of these material 
problems in HTSC cuprates are understood and as a re- 
sult consistent picture of tunnelling features is starting 
to emerge. 

From tunnelling experiments one obtains the energy 
gap in the superconducting state. Since we have already 
discussed this problem in [2], we will only briefly men- 
tion some important result, that in most cases, Gats' (V) 
has V-shape in all families of HTSC hole and electron 
doped cuprates. This is due to d — wave pairing with 
the gapless spectrum, which is definitely confirmed in 
the interference experiments on hole and electron doped 
cuprates [41]. Some experiments give the U-shape of 
Gns(V) which resembles s-wave pairing. This contro- 
versy is explained to be the property of the tunneling 
matrix element which filters out states with the maximal 
gap. 

Here we are interested in the electron-boson spectral 
function a 2 F(uo) in HTSC cuprates which can be ex- 
tracted by using tunnelling spectroscopy. We inform the 
reader in advance, that the shape and the energy width 
of a 2 F(u), which are extracted from the second deriva- 
tive d 2 I/dV 2 at voltages above the superconducting gap, 
in most HTSC cuprates resembles the phonon density 
of states F(u). This result is strong evidence for the 
importance of EPI in the pairing potential of HTSC 
cuprates. For instance, plenty of break-junctions made 
from Bi2212 single crystals [23| show that the negative 
peaks in d 2 I/dV 2 coincide with the peaks in the general- 
ized phonon density of states F p h(uj) measured by neu- 
tron scattering - see Fig. [20l 

The tunnelling spectra in Bi2212 break junctions (23[, 
which are shown in Fig. [20l indicates that the spectral 
function a 2 F(uj) is unchanged in magnetic field which 
disfavors SFI since in the latter case, this function should 
be sensitive to the magnetic field. The reported broaden- 
ing of the peaks in a 2 F(uo) are partly due to the gapless 
spectrum of d — wave pairing in HTSC cuprates. Addi- 
tionally, the tunnelling density of states Nt(ou) at very 
low T show a pronounced gap structure and it was found 
that 2A/T C = 6.2 - 6.5, where T c = 74 - 85 K and A is 
some average value of the gap. In order to obtain a 2 F{uo) 
the inverse procedure was<used by assuming s — wave 
superconductivity and the effective Coulomb parameter 
/i* « 0.1 [23]. The obtained a 2 F(uo) gives large EPI 
coupling constant X ep « 2.3. Although this analysis 
[23( ] was done in terms of s — wave pairing, it mimics 
qualitatively the case of d — wave pairing, since one ex- 
pects that d — wave pairing does not change significantly 
the global structure of d 2 I/dV 2 at eV > A albeit in- 
troducing a broadening in it - see the physical meaning 
in Appendix A. We point out that the results obtained 
in [23[ were reproducible on more than 30 junctions. In 
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FIG. 20: (a) Second derivative of I(V) for a Bi2212 break 
junction in various magnetic fields. The structure of minima 
in d 2 I/dV 2 can be compared with the phonon density of states 
F(oo); (b) the spectral functions a 2 F(oo) in various magnetic 
fields. From El 



that respect very important results on slightly overdoped 
Bi2212 — GaAs and on Bi2212 — Au planar tunnelling 
junctions are obtained in [24[ - see Fig. [2TJ 

These results show very similar features to those ob- 
tained in [23] on break-junctions. It is worth mentioning 
that several groups [25[, [26[, [2?J have obtained simi- 
lar results for the shape of the spectral function a 2 F{uo) 
from the I — V measurements on various HTSC cuprates 
- see the comparison in Fig. [22j The latter results leave 
no much doubts about the importance of the EPI in 
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FIG. 21: The spectral functions a 2 F(co) and the calculated 
density of states at OK (upper solid line) obtained from the 
conductance measurements the Bi(2212) — Au planar junc- 
tions. From [23] 



pairing mechanism of HTSC cuprates. 

In that respect tunnelling measurements on slightly 
overdoped Bi2Sr2CaCu20g [24|, [28] are impressive, 
since the Eliashberg spectral function a 2 F{uo) was ex- 
tracted from the measurements of d 2 I/dV 2 and by solv- 
ing the inverse problem. The extracted a 2 F(u) has sev- 
eral peaks in broad energy region up to 80 meV as it is 
seen in Fig. [2Ti \22\ which coincide rather well with the 
peaks in the phonon density of states F p h{oo). In [28[ nu- 
merous peaks, from PI — P13, in a 2 F(uo) are discerned 
as shown in FigJ23J which correspond to various groups 
of phonon modes - laying in (and around) these peaks. 
Moreover, in [24|, [28[ are extracted the coupling con- 
stants for these modes and their contribution (AT C ) to 
T c as it is seen in Fig. [23j Note, due to the nonlinearity 
of the problem, the sum of AT C is not equal to T c . 

The next remarkable result is that the extracted 
EPI coupling constant is very large, i.e. X ep (= 
2 / duoa 2 F(uo)/uo) = £\ A; « 3.5 - see Fig. [23j It is obvi- 
ous from Figs. (|23H24p that almost all phonon modes con- 
tribute to A ep and T c , which means that on the average, 
each particular phonon mode is not too strongly coupled 
to electrons thus keeping the lattice stable. None of the 
modes has too large A^ < 1.3, thus keeping the lattice 
stability. 

For a better understanding of the the EPI coupling in 
these systems we show in Fig. [25] the total and partial 
density of phononic states. In Fig.[23]it is seen that lower 
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FIG. 22: The spectral functions a 2 F(uj) from measurements 
of various groups: Vedeneev et al. 23f, G onnelli et al.[27j, 
Miyakawa et al. [25|, Shimada et al. [24|]. The generalized 
density of states GPDS for Bi2212 is at the bottom. From 



frequency modes from PI — P3, corresponding to Cu, Sr 
and Ca vibrations, are rather strongly coupled to elec- 
trons (with X K ~ 1) and give appreciable contributions to 
T c . It is also seen in Fig. [23] that the coupling constants 
Xi of the high-energy phonons (P9 — P13 with uo > 70 
meV) have A^ <C 1 and give moderate contribution to T c 
- around 10 %. These results confirm the importance of 
modes which cause the change of the Madelung energy in 
the ionic-metallic structure of HTSC cuprates, the idea 
also conveyed in [l| , Q - see more in Part II If definitely 
confirmed, these results are in accordance with the mod- 
erate oxygen isotope effect in cuprates near the optimal 
doping. We stress that each peak PI — P13 in a 2 F{uo) 
corresponds to many modes. In order to get filling on 
the structure of vibrations possibly strongly involved in 
pairing, we show in Figs. [26ll27l the structure of these 
vibrations at special points in the Brillouin zone. It is 
seen in Fig. [26] that the low- frequency phonons PI — P2 
are dominated by Cu, Sr, Ca vibrations. 

It is seen in Fig. [23] that P3 modes are stronger cou- 
pled to electrons than P4 ones, although the density of 
state for the P4 modes is larger. The reason for such 
an anomalous behavior might be due to symmetries of 
corresponding phonons as it is shown in Fig. [23 Namely 
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FIG. 23: Table I - Phonon frequency cj, EPI coupling constant 
Xi of the peaks PI — P13 and contribution AT C to T c of each 
peak in a 2 P(u;)-shown in Fig |24l obtained from the tunnelling 
conductance of Bi2Sr2CaCu20s- AT C is the decrease in T c 
when the peak in a 2 F(oo) is eliminated. From |28j ]. 



FIG. 24: The spectral functions a 2 F{uo) from the tunnelling 
conductance of P^SV^CaCi^Os for the positive and the neg- 
ative bias voltages, and the averaged one |24|. The averaged 
one is divided into 13 components. The origin of the ordinate 
is 2, 1, and -0.5 from the top down. From @, 



to the P3 peak contribute axial vibrations of 0(1) in 
the Cu2 plane which are odd under inversion, while in 
the P4 peak these modes are even. The in-plane modes 
of Ca and 0(1) are present in P3 which are in-phase 
and out-of-phase modes, while in P4 they are all out- 
of-phase modes. For more information on other modes 
P5 — P13 see [28]. We stress that the Eliashberg equa- 
tions based on the extracted a 2 F(uj) of the slightly over- 
doped Bi2Sr2CaCu20s with the ratio (2A/T C ) « 6.5 de- 
scribe rather well numerous optical, transport and ther- 
modynamic properties [28]. However, in underdoped sys- 
tems with (2A/T C ) « 10, where the pseudogap phenom- 
ena are pronounced, there are serious disagreements be- 
tween experiments and the Eliashberg-like theory. 

Similar conclusion, regarding the properties of the 
EPI spectral function a 2 F(uo) in HTSC cuprates, comes 
out from tunneling measurements on the Andreev (Z <C 
1, low barrier)- and Giaver (Z ^> 1, high barrier) -type 
jun ctions in La2- x Sr x CuO^ and YBCO compounds 
[971 ], where the extracted a 2 F(co) is in accordance with 
the phonon density of states F p h(co) see Fig. [28l 

Note that the BTK parameter Z is related to the trans- 
mission and reflection coefficients for the normal metal 
(1 + Z 2 )- 1 and Z 2 (l + Z 2 )' 1 respectively. 

Although most of the peaks in a 2 F(u) in HTSC 
cuprates coincide with the peaks in the phonon density of 
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FIG. 25: The phonon density of states F(u) (PDOS) of 
B%2Sr2C aCu20% compared with the generalized density of 
states (GPDOS) Atomic vibrations: 01 - O in the Cu0 2 
plane; 02 - apical O; 03 - O in the BiO plane. From [3 
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FIG. 26: Atomic polarization vectors and their frequencies (in 
meV) at special points in the Brillouin zone. The larger closed 
circles in the lattice are O-ions. r — X is along the Cu-O- 
Cu direction. Arrows indicate displacements. The modes in 
square and round brackets are the transverse and longitudinal 
optical modes respectively. (Top) - modes of the PI peak. 
(Bottom) - modes of the P2 peak. From [3, [HI. 



FIG. 27: Atomic polarization vectors and their frequencies 
(in meV) at special points in the Brillouin zone. The larger 
closed circles in the lattice are O-ions. T — X is along the Cu-O- 
Cu direction. Arrows indicate displacements. The modes in 
square and round brackets are the transverse and longitudinal 
optical modes, respectively. (Top) - modes of the P3 peak. 
(Bottom) - modes of the P4 peak. From [3, [HI. 



states it is legitimate to put the question - can the mag- 
netic resonance in the superconducting state give contri- 
bution to the a 2 F{uo)l In that respect very important 
inelastic neutron scattering measurements of the mag- 
netic resonance as a function of doping [98] give that the 
resonance energy E r scales with T c , i.e. E r = (5 — 6)T C 
as shown in Figl29l This means that if one of the peaks 
in a 2 F{uo) is due to the magnetic resonance at uo = E r , 
then it shifts strongly with doping as it is observed in 
(98[. This is contrary to phonon peaks (energies) whose 
positions are doping independent. To this end, recent 
tunnelling experiments on Bi2212 [29[ show clear doping 
independence of a 2 F(uo) as it is seen in Fig. [30l This 
remarkable result is an additional and strong evidence in 
favor of EPI and against the SFI mechanism of pairing 
in HTSC cuprates. 



It is interesting that in the vacuum tunneling STM 
measurements [99] the fine structure in d 2 I/dV 2 at eV > 
A was not seen below T c , while the pseudogap structure is 
observed at temperatures near and above T c . This result 
could mean that the STM tunnelling is likely dominated 
by the nontrivial structure of the tunnelling matrix ele- 
ment (along the c-axis), which is derived from the band 
structure calculations JlOOll . However, recent STM ex- 
periments on Bi2212 [30[ give important information on 
possible nature of the bosonic mode which couples with 
electrons. In [30] the local conductance dI/dV(r,E) is 
measured where it is found that d 2 I/dV 2 (r,E) has peak 
at E(r) = A(r) +ft(r) where dI/dV{r,E) has the maxi- 
mal slope - see Fig. Icffla). 

It turns out that the average phonon energy Cl depends 

- —1/2 

on the oxygen mass, i.e. ft ~ M Q 1 , with Qiq = 52 meV 
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FIG. 28: (a) d 2 I/dV 2 of a Giaver-like contact in La 
2-xStxCuOa, - note the large structure below 50rael/; (b) 
d 2 I/dV 2 of an Andreev- and Giaver-like contact compared to 
the peaks in the phonon density of states. From [97j] 



and £li8 ~ 48 meV - as it is seen in Fig. l3lfb). This re- 
sult is a convincing evidence that phonons are strongly 
involved in the quasi-particle scattering. A possible ex- 
planation is put forward in [30[ by assuming that this 
isotope effect is due to the phonon which interacts 
with anti- nodal quasi-particles. 



In our opinion the important message of tun- 
nelling experiment s in HTSC cuprates (by including 
Bai- x K x BiOs too |lOlj , |lQ2| ) is that there is strong ev- 
idence for the importance of EPI and that no particular 
phonon mode can be singled out in the spectral function 
a 2 F{uj) as being the only one which dominates in pair- 
ing mechanism. This important result means that the 
high T c is not attributable to a particular phonon mode 
in the EPI mechanism, i.e. all phonon modes contribute 
to X ep . Having in mind that the phonon spectrum in 
HTSC cuprates is very broad (up to 80 meV ), then 
the large EPI constant (X ep > 2) obtained in tunnelling 
experiments is not surprising at all. 



140 : 
120 ; 
100 - 

80 7 

60 P 

40 

20 h 



h — i — i — i — |~~n — m — | — n — r i [ t i i t _ 
— O-YBCO-E 



E r (meV) 



BSCO - E 



(SIS) 



5,3 k T 

,B c 



Underdoped Overdoped 

J— l I I I I k 1 L I I I t I I 



150 



100 H 

o 

50 w 



-0.1 -0,05 



0.05 0.1 



n -n 

h opt 



FIG. 29: Doping dependence of the energy E r of the magnetic 
resonance peak at ir, tt in YBCO and Bi2212 measured at low 
temperatures by inelastic neutron scattering. From [98[ 
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FIG. 30: Second derivative of I(V) for a Bi22l2 tunnelling 
junctions for various doping: UD-underdoped; OD-optimally 
doped; OVD-overdoped system. The structure of minima in 
d 2 I/dV 2 can be compared with the phonon density of states 
F(ui). The full and vertical lines mark the positions of the 
magnetic resonance energy E r ~ 5.4T C for various doping 
taken from Fig |29l Red tiny arrows mark positions of the 
magnetic resonance E r in various doped systems. Dotted ver- 
tical lines mark various phonon modes. From [29[ 
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FIG. 31: (a) Typical conductance dl / dV (r , E) . The ubiq- 
uitous feature at eV > A (gap) with maximal slopes, which 
give peaks in d 2 I/dV 2 (r, E) are indicated by arrows, (b) The 
histograms of all values of H(r) for samples with 1 6 - right 
curve and with 1 8 - left curve. From [30] 



E. Phonon spectra and EPI 

Although experiments, such as inelastic neutron and 
Raman scattering, related to phonon spectra and their 
renormalization by EPI do not give directly a 2 F(cj), 
as the tunnelling and optic spectra do, they neverthe- 
less give useful information on the strength of EPI for 
some particular phonons. We stress in advance that the 
interpretation of experimental results in terms of EPI 
for weakly interacting electrons might be risky since in 
strongly correlated systems, such as HTSC cuprates, the 
phonon renormalization due to EPI is rather different 
from that in weakly correlated metals - see [103j . Since 
these questions are thoroughly studied in the excellent 
review |lQ3| - see also Part II, we shall briefly enu- 
merate the main points: (1) In strongly correlated sys- 
tems, the EPI coupling for a number of phononic modes 
can be significantly larger than the LDA and Hartree- 
Fock met hods predict. This is due to many-body ef- 
fects |lQ3( . (2) In strongly correlated systems the quasi- 
particle charge susceptibility, which enters the phonon 
self-energy II (q, u)(~ x c (q, co ))^ * s milcn more suppressed 
than in weakly correlated metals and these effects are out 
of the LDA possibilities 0, [lQ3|, [lQjj. This is one of the 
reasons that the analysis of experiments on phonon renor- 
malizations in the framework of LDA underestimates the 
EPI coupling constant significantly - on all these ques- 
tions, see more details in Part II 



1. The phonon Raman scattering 

The phonon Raman scatteri ng gi ve s als o e vide nce for 
appreciable EPI in cuprates |lll| ], |l!2| , Ill3( . We 
enumerate some of them - see more in Q and Ref- 
erences therein, (i) There is a pronounced asymmet- 



ric line-shape (of the Fano resonance) in the metallic 
state. For instance, in Y Ba^Cu^O^ two Raman modes 
at 115 era -1 (Ba dominated mode) and at 340 cm -1 (O 
dominated mode in the Cu02 planes) show pronounced 
asymmetry which is absent in YBa2CusOe. This re- 
sult points to appreciable interaction of Raman activ e 
phon ons with continuum states (quasi-particles) jlll| . 
[U2j,. (ii) The phonon frequencies for some A\ g and B\ g 
are strongly renormalized in the superconducting state, 
between (6 — 10) %, pointing again to large EPI [l!3| ] 
- see also in 0, [44j |. To this point we mention that 
the electronic Raman scattering in cuprates show a re- 
markable correlation between the Raman cross-section 
^exp(^) and the optical conductivity in the a-b plane 
ct(cj), i.e. S exp (uj) ~ a(uo) [2]. It was argued above 
that EPI with the very broad spectral function a 2 F{uo) 
explains in a natural way the uo and T dependence of 
<j(lj). This means that the electronic Raman spectra in 
cuprates can be explained by EPI in conjunction with 
strong correlations. This conclusion is su pported by cal- 
culations of the Raman cross-section |l05| which take into 
account EPI with a 2 F{uo) extracted from tunnelling mea- 
surements in YBa2Cus06-\- x and B%2Sr2CaCu20s+x 0, 
[23[ -J28[. Quite similar properties (to cuprates) of the 
electronic Raman scattering, as well as of a(w), R(uo) 
and p(T), were observed in experiments [(0 on isotropic 
3D metallic oxides Lao.^Sro^CoOs and Cao^Sro^RuOs 
where there are no signs of antiferromagnetic fluctua- 
tions. This means that low-dimensionality and antifer- 
romagnetic spin fluctuations cannot be a prerequisite for 
anomalous scattering of quasi-particles and EPI must be 
inevitably taken into account since it is present in all 
these compounds. 



2. Neutron scattering, phonon spectra and EPI 

The softening of numerous phonon modes has been ob- 
served in the normal state of cuprates giving important 
evidence for pronounced EPI. There are several impor- 
tant reviews on this subject [108] and here we discuss 
briefly two important examples which demonstrate in an 
impressive way the inefficiency of the LDA band struc- 
ture calculations to treat quantitatively and qualitatively 
EPI in HTSC cuprates. Namely, by doping the Cu-0 
bond-stretching phonon mode shows a substantial soften- 
ing at q^6 = (0.5, 0, 0) - called the half-breathing phonon, 
and a large broadening by 5 meV at 15% doping. While 
the s oftening can be partly described by the LDA method 
[109|, it predicts an order of magnitude smaller broad- 
ening than the experimental one. The reason for this 
failure lies in strong correlations, which are not in c luded 
in the LDA method as explained in 0, @, [ToJ, [lO^ . 
They give rise to an increase of the EPI coupling and to 
strong suppression of the charge fluctuations which en- 
ter the phonon self-energy via the charge susceptibility 
- see below and in Part II The neutron scattering in 
Lai, 85 Sro.isCuO^ gives evidence for large (30%) soften- 
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ing of the Of with Ai symmetry with the energy uo « 60 
meV, which is theoretically predicted in [106], and for 
the large line- width about 17 meV which also suggest 
strong EPI. 

As it is discussed in the Introduction, there are re- 
cently several calculations of the EPI coupling constant 
X ep in the framework of DFT (or LDA), where very 
small X ep « 0.3 was obtained [l3j], (TJ. However, the 
LDA method is inadequate for strongly correlated sys- 
tems as it does not correctly take into account exchange- 
correlations and many-body effects, and therefore over- 
estimates the screening in cuprates. If DFT is able to 
describe EPI correctly, it must also be able to calculate 
phonon renormalization, such as softening and broaden- 
ing of the spectrum. In fact DFT completely fails to de- 
scribe this renormalization for some important phonon 
modes and therefore fails to describe the effect of EPI on 
the electronic spectrum. The critique of LDA (DFT) 
results in HTSC cup rate s is done in [2] and recently 
strongly argued in Q , |l!0( by its disagreement with neu- 
tron scattering measurements as it is shown in Fig. [32l 

The point is that DFT (LDA) can reproduce the 
phonon softening in Lai. 85 Sr 0.15CUO4 and Yb^Cu^Oj 
rather good at low momenta q = (h, 0, 0) but predicts 
smooth softening at higher q, while experiments show 
pronounced features for h = 0.3. At the same time DFT 
predicts an order of magnitude smaller line width than 
experiments |l!4j . 

In Part II we shall discuss some theoretical approaches 
related to EPI renormalization of phonons in strongly 
correlated systems. Here, we point out two results. 
First, there is an appreciable difference in the phonon 
renormalization in strongly and weakly correlated sys- 
tems. Namely, the change of phonon frequencies in the 
presence of conduction electrons is proportional to the 
coupling constant |g q | and charge susceptibility x c > i- e - 



Mq) ~ Itfql Re X c 

l#q| 2 l Im Xcl 



«>(q) 



while the line width is given by 
It turns out that in strongly cor- 
related systems doped with hole concentration ^ « 1 
the charge fluctuations are s uppr essed in which case the 
following sum-rule holds |l04j 



ttN ^ / 



\lm X M\du; = 2S(l-S)N, 



70 
65 



T~i — 1 — 1—1 — 1 — 1 — 1— 1 — 1 — r 

La, K Si 01% CuO 4 




-OFT (Giusinoet al y 45 



YBa CuA 




0.0 1 2 03 0.4 05 
h(r. Iai.) 



q=iOhO> 



-DFT <Bohranet al.) 



0.0 0.1 0.2 3 0.4 05 

h(r.l.u.) 



16 
14 
12 

§ B 

i 6 
4 

2 





• *. 

■„•• 

• DFT {Bohnenot ah) 
m Ex penny nl ^ 







0.2 4 
h(r.l.u.) 



FIG. 32: Comparison of DFT calculations with experimen- 
tal results:(a) in Lai.ssSro.isCuOA] (b) in YBaiCu^Oi . (c) 
Phonon line widths in La\.%*>Sr§.\hCuO±. DFT calculations 
1311 gives much smaller width than experiments [1141 ]. From 



not present in the latter |l04l ]. It can be shown that 



<l\LDA' 



for some phonon modes one has |g q | ^> \g t 
instance, for the half-breathing mode, one has |g q | 



For 



3 \gd\ LDA that is first calculated in [104[. These results 
point to inadequacy of LDA in calculations of EPI effects 
in HTSC cuprates. 



while in the LDA method one has 

\LDA 



1 poo 



du = (l- 5)N. 



This means that for low doping J C 1 (note n = 1 — 5), 
one has |Imx c | <C \^- m Xc DA \ an d LDA strongly under- 
estimates the coupling constant, i.e. \g<i\ LDA <^ |g q |. We 
stress that there is no such strong suppression in the 
quasi-particle self-energy [5[. 

Second, the theory gives that the coupling constant 
|<7 q | in HTSC cuprates can be significantly larger than 
LDA predicts, which is due to some many-body effects 



F. Isotope effect for various doping 

The isotope effect olt c in the critical temperature T c 
was one of the very important proof for the EPI pairing in 
low-temperature superconductors (LTSC). As a curiosity 
the isotope effect in LTSC systems was measured almost 
exclusively in monoatomic systems and in few polyatomic 
systems: the hydrogen isotope effect in PdH, the Mo 
and Se isotope shift of T c in MoqSc%, and the isotope ef- 
fect in NbsSn and MgL>2. We point out that very small 
(ctT c ~ in Zr and Ru) and even negative (in PdH) 
isotope effect in some polyatomic systems of LTSC ma- 



30 



terials are compatible with the EPI pairing mechanism 
but in the presence of substantial Coulomb interaction or 
lattice anharmonicity. The isotope effect ar c cannot be 
considered as the smoking gun since it is sensitive to nu- 
merous influences. For instance, in Mgl>2 it is with cer- 
tainty proved that the pairing is due to EPI and strongly 
dominated by the boron vibrations, but the boron isotope 
effect is significantly reduced, i.e. ar c ~ 0.3. It is still 
unexplained. The situation in HTSC cuprates is much 
more complicated because they contain many-atoms in 
unit cell. Additionally, the situation is complicated with 
the presence of intrinsic and extrinsic inhomogeneities 
which can mask real effects. On the other hand new 
techniques such as ARPES, STM, jiSR allow studies of 
the isotope effects in quasi-particle self-energies, i.e. as, 
which will be discussed below. 



x = 0.125 with T c <C T c>max one has ac u ~ 0.8 — 1 
with olcu + olo ~ 1.8 [107j. The appreciation of copper 
isotope effect in La2- x M x Cu04 tells us that vibrations 
other than oxygen ions are important in giving high T c . 
In that sense one should have in mind the tunnelling ex- 
periments discussed above, which tell us that all phonons 
contribute to the Eliashberg pairing function a 2 F(k, uj) 
and according to these results the oxygen modes give 
moderate contribution to T c [28[. Having these facts in 

mind, then the small oxygen isotope effect a^ in opti- 
mally doped cuprates, if it is intrinsic property, does not 
exclude the EPI mechanism of pairing. 



2. Isotope effect as in the self-energy 



1. Isotope effect ar c in T c 

This problem will be discussed only briefly since more 
extensive discussion can be found in [2] . It is well known 
that in the pure EPI pairing mechanism, the total isotope 
coefficient a is given by 

a T =V Q M = -V dlnTc (41) 

l,p l,p I 

(p) 

where M\ is the mass of the i-th element in the p-th 
crystallographic position. Note that in the case when the 
screened Coulomb interaction is negligible, i.e. //* = 0, 
one has olt c — 1/2- From this formula one can de- 
duce that the relative change of T c , 5T C /T C1 for heav- 
ier elements is rather small - for instance it is 0.02 
for 135 Ba -^ 138 Ba, 0.03 for 63 Cu -^ 65 Cu and 0.07 
for 138 La -^ 139 La. This means that measurements 
of OLi for heavier elements are at /or beyond the abil- 
ity of present day experimental techniques. Therefore 
most isotope effect measurements were done by substi- 
tuting light atoms 16 O by 18 O only. It turns out that 
in most optimally doped HTSC cuprates ao is small. 
For instance ao ~ 0.02 — 0.05 in YBCL2CU3O7 with 
T C)max « 91 K, but it is appreciable in La1.s55ro.15Ct/O4 
with T C5max ~ 35 K where ao ~ 0.1 — 0.2. In 
Bi2Sr2CaCu20s with T CjYnax w 76 K one has ao ~ 
0.03 — 0.05 while ao ~ 0.03 and even negative (—0.013) 
in £/25r 2 Ca 2 C?/20io with T c , max « 110 K. The ex- 
periments on Tl2Ca n -\BaCu n 02n+A = 2,3) with 
L C5m ax ~ 121 K are still unreliable and ao is unknown: 
In the electron-doped (Nd\- x Ce x )2CuO^ with T C)ITiax « 
24 K one has ao < 0.05 while in the underdoped materi- 
als ao increases. The largest ao is obtained even in the 
optimally doped compounds like in systems with substi- 
tution, such as Lai.855ro.i5Ci£i_ x M £C 04, M = Fe,Co, 
where ao ~ 1.3 for x « 0.4 %. In La2- x M x Cu04 there 
is a Cw isotope effect which is of the order of the oxy- 
gen one, i.e. ac u ~ <^o giving ac u + <^o ~ 0.25 — 0.35 
for optimally doped systems (x = 0.15). In case when 



The fine structure of the quasi-particle self-energy 
E(k, cj), such as kinks, can be resolved in ARPES mea- 
surements and in some respect in STM. It turns out 
that there is an isotope effect in the se lf-en er gy in the 
optimally doped 5/2212 samples [H|, EH, VlM- In 
the first paper on this subject |8l| , there is a red shift 
b~Uk,7o ~ —(10 — 15) meV of the nodal kink at LJkjo — 70 
meV for the 16 O ^ 18 O substitution. This isotope shift 
of the self-energy ST, = T>iq — T>is ~ 10 meV is more 
pronounced at large energies uo = 100 — 300 meV . How- 
ever, there is a dispute on the la tter result which is not 
confirmed experim ental ly [l 15ll , |l!6j ]. The isotope ef- 
fect in ReS(k,cj) [HI, II 16 11 can be well described in 
the framework of the Migdal-Eliashberg theory for EPI 
[82| which is in accordance with the recent ARPES mea- 
surements with low-energy photons ~ 7 eV jll7j ]. The 
latter allowed very good precision in measuring the iso- 
tope effect in the nodalpoint of Bi2212 with T c 16 = 92.1 
K and T c 18 = 91.1 K [mj- They observed shift in the 
maximum of ReE(k7v, uo) - at 0^,70 w 70 meV which cor- 
responds to the half-breathing or breathing phonon, by 
b~Uk,7o ~ 3.4 ± 0.5 meV as shown in Fig. [33j 

By analyzing the shift in ImX!(kA/-, uj) - shown in Fig. 
[33l one finds similar result 6uk t 70 ~ 3.2 ± 0.6 meV. The 
similar shift was obtained in STM measurements [3O] 
which is shown in Fig. l3Tlfb) and can have its origin in 
different phonons. We would like to stress two points: (i) 
in compounds with T c ~ 100 K the isotope effect in T c 
is moderate, i.e. aL ' < o.i, 

If we consider this 

value to be intrinsic then it is not in conflict with the 
tunnelling experiments in [28] which give evidence that 
vibrations of heavier ions contribute significantly to T c - 
see the above discussion on the tunnelling spectroscopy; 
(ii) the extracted value in |ll7j of the effective EPI cou- 
pling constant ~ 0.6, which is smaller than the real 
X ep - see above the discussion on ARPES, is significantly 
larger than the LDA theory predicts X e ^ DA < 0.3 [l3[, 
jl4j ]. This again points that the LDA method does not 
pick up the many-body effects due to strong correlations- 
see Part II. 
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FIG. 33: (a) Effective ReY> for five samples for O 16 (blue) and 
O 18 (red) along the nodal direction, (b) Effective ImS deter- 
mined from MDC full widths. An impurity term is subtracted 
at uj = 0. From [HI- 



IV. SUMMARY OF PART I 

The analysis of experimental data in HTSC cuprates 
which are related to optics, tunnelling and ARPES mea- 
surements near the optimal doping give evidence for the 
large electron-phonon interaction with the coupling con- 
stant 1 < X ep < 3.5. The analysis is done in the frame- 
work of the Migdal-Eliashberg theory which is a reliable 
approach for systems near the optimal doping. The spec- 
tral function averaged over the Fermi surface a 2 F{uo) is 
extracted from various tunnelling measurements and it 
contains peaks at the same positions as the phonon den- 
sity of states. The energy width of a 2 F(u) is the same 
as the width of the phonon density of states F(w). This 
is an unambiguous proof for the important role which 
EPI plays in the pairing mechanism of cuprates. The op- 
tical IR reflectivity data provide additional support for 
this finding since the transport spectral function has the 
width and global properties similar to F(ou). These find- 
ings are additionally and strongly supported by ARPES 
measurements on BISCO compounds. The ARPES kinks 
in the quasi-particle self-energy can be explained exclu- 
sively by EPI and there is no much room for spin fluctu- 
ations (SFI) effects. The weakness of SFI is unambigu- 



ously proved in magnetic neutron scattering on YBCO 
where the imaginary part of the susceptibility is dras- 
tically reduced in the low energy region by going from 
slightly under doped toward optimally doped systems, 
while T c is practically unchanged. This implies that the 
SFI coupling constant is limited to the value X s ^ < 0.3. 
All these results do not leave doubts on the significance 
of EPI and weakness of SFI. The obtained total EPI cou- 
pling constant is rather large, i.e. 1 < X ep < 3.5, while 
the transport coupling constant is X tr ~ A/3. The differ- 
ent renormalization of the quasi-particle and transport 
self-energies by the Coulomb interaction and strong cor- 
relations points to dominance of the small-momentum 
scattering in EPI. This will be discussed in Part II 

Inelastic neutron scattering measurements in cuprates 
show that the broadening of phonon lines is by an order 
of magnitude larger than the LDA (DFA) method pre- 
dicts. Since the phonon line-widths depend on the EPI 
coupling and the charge susceptibility it is evident that 
calculations of both quantities are beyond the range of 
applicability of LDA. As a consequence, LDA overesti- 
mates electronic screening and thus underestimates the 
EPI coupling. This means that LDA is suitable only for 
weakly correlated systems, while many-body effects due 
to strong correlations are not contained in this mean-field 
type theory. 

In spite of the very promising and encouraging results 
about the dominance of EPI in cuprates the theory is still 
confronted with the task of obtaining sufficiently large 
coupling constant in the d-channel in order that EPI con- 
forms with d-wave pairing. At present we do not have 
such a detailed microscopic theory although some con- 
cepts, such as the the dominant EPI scattering at small 
transfer momenta, are understood at least qualitatively. 
These set of problems and questions will be discussed in 
Part II 
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V. APPENDIX - SPECTRAL FUNCTIONS 

A. Spectral functions a 2 F(k, k', uj) and a 2 F{uj) 

Before discussing experimental results in cuprate su- 
perconductors which are related to various spectral func- 
tions, we introduce the reader briefly into the sub- 
ject. The quasiparticle-bosonic (Eliashberg) spec- 
tral function a 2 F(k, k 7 , uj) and its Fermi surface aver- 
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age a 2 F(u) = ((a 2 F(k, k', ^))) k k , enter the quasi- 
particle self-energy E(k, a;), while the transport spec- 
tral function a 2 F tr (uj) enters the transport self-energy 
Dtr(k, oj) and dynamical conductivity <j(oj). Since 
the Migdal-Eliashberg theory for EPI is well defined 
we define the spectral functions for this case and 
the generalization to other electron-boson interaction 
is straightforward. In the superconducting state the 
Matsubara Green's function G(k, uj n ) and E(k, uo n ) are 
2x2 matrices with the diagonal elements G\\ = 
G(k, u; n ), G22, ^ii = E(k, (j n ), S22 and off-diagonal el- 
ements G12 = F(k,^ n ),G 2 i,Hi2 = $(k,o; n ),E2i 
which describe superconducting pairing. By defin- 
ing [1 - Z(k, cj n )] = [E(k, uj n ) - E(k, -w n )] /2 and 
x(k,cj n ) = [E(k,a; n ) + E(k, — o; n )]/2, the Eliashberg 
functions for EPI in the presence of the C oulo mb inter- 
action (in the pairing channel) read |l!8| . |ll9^ . |l20j 



Z(k 



'" nJ N{ji)uj n D(p,u, m y 



p,m 



(42) 



X(k,^ n ) 



$(k,cj n ) 



Ag p (^ n - o; m ) e (p) - ji + x(p, cj m ) 



iV(/i) 



£>(p,^m) 



(43) 



N ^ 

p,m 



A kp 







AT(/i) 



k P 



$(p,^m) 



Z>(p,cj m )' 
(44) 

where iV(/i) is the density of states at the Fermi sur- 
face, uo n = 7rT(2n + l), <I>(k, a; n ) = Z(k, c^ n )A(k, c^ n ) and 
L> = cj^Z 2 + (e — fi + x) 2 + $ 2 - (For studies of optical 
properties - see below, it is useful to introduce the renor- 
malized frequency iuj n (iuj n ){= iuj n Z(uj n )) = uj n — E(o; n ) 
(or its analytical continuation uj(uj) = Z(uj)uj = uo — 
E(o;)). These equations are supplemented with the elec- 
tron number equation n(/i) (/i is the chemical potential) 

2T 



™0) = — ^G(p,u; m )e* 



where B K (k — p;i>) is the phonon spectral function of the 
ft-th phonon mode related to the phonon propagator 



f 

Jo 



-B K (w)dv. (48) 



The renormalized coupling constant #J7kp(~ 9k kp £_1 7) 
comprises the screening effect due to long-range Coulomb 
interaction (~ e~ x - the inverse electronic dielectric func- 
tion) and short-range strong correlations (~ 7 - the ver- 
tex function) - see more in Part II Usually in the case 
of low-temperature superconductors (LTS) with s-wave 
pairing the anisotropy is rather small (or in the presence 
of impurities it is averaged out) whi ch allo w s an averag- 
ing of the Eliashberg equations [ll8( ] , |l!9( , |l2(| 



Z{uo n ) 



X(LU r , 



A 2 (w m )' 



(49) 



Z(w„)A(w„) =ttT^[A( 



- n(uO c )0(uO c ~ \uO r 



A(o; m ) 



A(c^ n - U0 m ) 



f 

Jo 



dv 



V^ + A 2 (o; m ) 

2^a 2 F(^) 
i/ 2 + (cj n - CJ m ) 2 ' 



(50) 



(51) 



Here a 2 F(c^) = ((a 2 F(k, k', w))) k>k/ where ((...)> k>k / is 
the average over the Fermi surface. The above equations 
can be written on the real axis by the analytical contin- 
uation iuo m — > uo + iS where the gap function is complex 
i.e. A(uj) = Ar(oj) -\-iAj(uj). The solution for A (a;) al- 
lows the calculation of the current- voltage characteristic 
/(V) and tunnelling conductance Gns(V) = dlws/dV in 
the superconducting state of the NIS tunneling junction 
where Ins(Y) * s gi yen by 



I NS (V) = 2eJ2\T Kp \ 2 I 
k, P J ~ 



™ du 
-00 27r 



1 - — e ( p ) ~^ + *( p ^) (45) 



Note that in the case of EPI one has A kp (^ n ) = 

^k P (^n)(= Akp(^n)) (with v n = irTn) where A kp (^ n ) 
is defined by 



Akp(^n) = 2 



2/2 + "I 



(46) 



Ajv(k,o;)A s (p,a; + e^)[/M - f(uo + eV)]. (52) 

Here, A/v^k, a;) = — 2/mG j /v,s , (k, a;) are the spectral 
functions of the normal metal and superconductor, re- 
spectively and f(uo) is the Fermi distribution function. 
Since the angular and energy dependence of the tun- 
nelling matrix elements | Tk, p | 2 is practically unim- 
portant for s — wave superconductors, then in that case 
the relative conductance ctns(V) = Gns(V)/Gnn(V) is 
proportional to the tunnelling density of states Nt(oj) = 
J A s (k,oj)d 3 k/(27r) 3 , i.e. (t ns (oj) « N t (uj) where 



al p F(v) = Nfjj) to P | B ^ k ~ P^) ( 47 ) 



A^t(cj) = i?e 



/ (cj + z7(cj)) 2 -Z 2 (^)A(cj) 2 



(53) 
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Here, Z(u) = Z(u)/ReZ(uj), = j(u)/ReZ(uj), 

Z(uj) = ReZ(uo) + ij(uj) /w and the quasi-particle scat- 
tering rate in the superconducting state 7 s (cj,T) = 
-2/mS(w,T) is given by 



7 s (cj,T) = 2tt / dva 2 F(v)N s (v + u){2n B (v) 



+ n F (y + a;) + n F (^ - cj)} + 7*' 



(54) 



where N s (us) = Re{o;/(a; 2 — A 2 ) 1 / 2 is the quasi-particle 
density of states in the superconducting state, ub,f(v) 
are Bose and Fermi distribution function respectively. 
Since the structure of phononic spectrum is contained 
in a 2 F(uj), it is reflected on A (a;) for uo > Ao (the real 
gap obtained from A = ReA(uo = Ao)) which gives the 
structure in Gs(V) at V = Aq+o;^. On the contrary one 
can extract the spectral function a 2 F(uo) from Gns (V) 
by the inver sion procedure proposed by McMillan and 
Rowell jl2l| . Q. It turns out that in low-temperature 
superconductors, negative peaks of d 2 I/dV 2 at eVi = 
A + 00ph.ii correspond to the peak positions of a 2 F(uj) 
and F{uj). However, we would like to point out that 
in HTSC cuprates the gap function is unconventional 
and very anisotropic, i.e. A(k, iuo n ) ~ cos k x a — cos k y a. 
Since in this case the extraction of a 2 F(k, k 7 , ui) is ex- 
tremely difficult and at present rather unrealistic task, 
then an "average" a 2 F{uo) is extracted from the experi- 
mental curve GsiV). There is belief that it gives relevant 
information on the real spectral function such as the en- 
ergy width of the bosonic spectrum (0 < uo < o; max ) and 
positions and distributions of peaks due to bosons. It 
turns out that even such an approximate procedure gives 
valuable information in HTSC cuprates - see discussion 
in Section III D. 

Note that in the case of spin-fluctuation interaction 
(the SFI model) one should make difference between 
Akp(i^n) and X^ p (iu n ) since they differ by sign i.e. 

Akp(^n) = — A^p(^n) > since SFI is repulsive in the 
pairing-channel - see Eqs. ([38]). 



1. Inversion of tunnelling data 

Phonon features in the conductance crNsiY) at eV = 
Ao + ujph makes the tunnelling spectroscopy a power- 
ful method in obtaining the Eliashberg spectral function 
a 2 F{uo). Two methods were used in the past for extract- 
ing a 2 F(u). 

The first method is based on solving the inverse prob- 
lem of the nonlinear Eliashberg equations. Namely, by 
measuring (Tns(V)i one obtains the tunnelling density of 
states Nt(oo)(~ (Tns(^)) and by the inversion procedure 
one gets a 2 F(uj) jl2l| . In reality the method is based on 
the iteration procedure - the McMillan-Rowell (MR) in- 
version, where in the first step an initial a 2 Fi n i(cu), \i\ ni 



and Ai n i(uj) are inserted into Eliashberg equations (for 
instance A ini (uj) = Ao for uo < uo and A^(c^) = for 
uo > ouq) and then (Ji n i(V) is calculated. In the second 
step the functional derivative 5a(uo) / 'da 2 F{uo) (uo = eV) 
is found in the presence of a small change of a 2 Fi n i(uo) 
and then the iterated solution a 2 F^(uo) = a 2 F ini (uj)-\- 
5a 2 F(u) is obtained, where the correction 5a 2 F(uo) is 
given by 



5a 2 F(uo) 



dv ^ J^p^ 1 [o~exp(v) ~ o- ini (is)}. (55) 



The procedure is iterated until a 2 F^{uo) and /ij^ con- 
verge to a 2 F{uo) and /i*which reproduce the experimen- 
tally obtained conductance cr e ^g(V). In such a way the 
obtained a 2 F{uo) for Pb resembles the phonon density 
of states Fpb(u), that is obtained from neutron scatter- 
ing measurements. Note that the method depends ex- 
plicitly on /i* but on the contary it requires only data 
on ctns(V) up to the voltage Vm ax = ^£}f x + Ao where 



uo. 



ph 



is the maximum phonon energy (a F(uu) = for 
uj > ^JJ)f x ) and Ao is the zero-temperature superconduct- 
ing gap. One pragmatical feature for the interpretation 
of tunnelling spectra (and for obtaining the spectral pair- 
ing function a 2 F(uo)) in LTS and HTSC cuprates is that 
the negative peaks of d 2 I/dV 2 are at the peak positions 
of a 2 F(oj) and F(uu). This feature will be discussed later 
on in relation with experimental situation in cup rates. 

The second method has been invented in |l22| and it 
is based on the combination of the Eliashberg equations 
and dispersion relations for the Greens functions - we call 
it GDS method. First, the tunnelling density of states 
is extracted f rom the tunneling conductance in a more 
rigorous way |l23^ 



N T (V) 



1 



0~NS(V) _ 

<tnn(V) <t*(V) 
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Jo 



du 



x [N T (V - 

du 



u) - N T (V)} 



(56) 



where cr*(V) = exp{— /3V}o-nn(V) and the cons tant (3 is 
obtained from a n n (V) at large biases - see [ 1 2 2| . N T (V) 
under the integral can be replaced by the BCS density 
of states. Since the second method is used in extracting 
a 2 F{uo) in a number of LTSC as well as in HTSC cuprates 
- see below, we describe it briefly for the case of isotropic 
EPI a t T= Q K. I n that case the Eliashberg equations are 

E3, nn, mi, mi 



Z(uj)A(uj) 
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J A, 



duj' Re 



A(a/) 



A 2 (co')} 



1/2 



x [K + (lo, to')} — /i*9(lu c — w) 



(57) 
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duj' Re 



[w' 2 - A 2 (w')] 
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K-{w,u') 
(58) 
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where 
K± 



dva 2 F(y)( - 



1 



uj' — co + v — iO+ 



)• 



(59) 



Here /i* is the Coulomb pseudopotential, the cutoff uj c is 
approximately (5 — 10) cj^ x , Ao = A(Aq) is the energy 
gap. Now by using the dispersion relati on fo r the matrix 
Greens functions G(k, uj n ) one obtains |l22| 



Im S(uj) 



2uj 



duo 



f N T (uJ f )-N B cs(^) 



(60) 



1/2 

where S(u) = uj/ [oj 2 — A 2 (u)] . From Eqs. 
one obtains 



poj — A q 

J dva 2 F(uo-v)Re^A(v) [v 2 - A 2 {y)\ 1/2 } 



x ; 1 dva 2 F{v)N T (uo - v) + — 



UJ 



ImA(cj) 



/ duj'N T {uj f ) 
Jo Jo 



J p h 7 2a 2 F(u) 
dv- 



Based on Eqs. 
tracting a 2 F(uj) 



(u' + ^) 2 - UJ 2 ' 
(61) 

one obtains the scheme for ex- 



VNs(y),* NN (V)^N T (y) 



-> Im5(w) -> A(w) a 2 F(cj). 

The advantage in this method is that the explicit knowl- 
edge of /i* is not required |l22| . However, the integral 
equation for a 2 F(uj) is linear Fredholm equation of the 
first kind which is ill-defined - see the duscussion in Sec- 
tion II.B.2. 



2. Phonon effects in Nt{w) 

We briefly discuss the physical origin for the phonon 
effects in Nt(oo) by considering a model with only one 
peak, at c^o, in the phonon density of states F(uj) by as- 
suming for simplicity /i* = and neglecting the weak 
structure in Nt(oo) at ncoo + Ao, which is du e to the 
nonlinear structure of the Eliashberg equations [125| . In 
Fig EH it is seen that the real gap Ar(uj) reaches a max- 
imum at c^o + Ao then decreases, becomes negative and 
zero, while A/ (lj) is peaked slightly beyond ujo + Ao that 
is the consequence of the effective electron-electron inter- 
action via phonons. 
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FIG. 34: (a) Model phonon density of states F(oo) with the 
peak at ujq. (b) The real (solid) Ar and imaginary part A/ 
(dashed) of the gap A(cj). (c) The normalized tunnelling 
density of states Nt(w)/N(0) (solid) compared with the BCS 
density of states (dashed). From [l25l ]. 



It follows that for uj < ujo + Ao most phonons have 
higher energies than the energy uj of electronic charge 
fluctuations and there is over-screening of this charge by 
ions giving rise to attraction. For uj « ujo + Ao charge 
fluctuations are in resonance with ion vibrations giving 
rise to the peak in Ar(u). For c^o + Ao < uo the ions move 
out of phase with respect to charge fluctuations giving 
rise to repulsion and negative Ar{uj). This is shown in 
Fig. 1531b). The structure in A (a;) is reflected on Nt(uj) 
as shown in Fig. l34lfc) which can be reconstructed from 
the approximate formula for Nt(oj) expanded in powers 
of A/uj 



N t (uj) 



1 

1+ 2 



A r (uj) 



A t (uj) 



this gives Nt(w) > 
the real value Ar(oj) 



As Ar(uj) increases above Ao 
Nbcs(u), while for uj > c^o + Ao 
decreases while Ai(co) rises and Nt(oo) decreases giving 
rise for Nt(uj) < Nbcs(u). 

B. Transport spectral function o% r F(ui) 

The spectral function a 2 r F{uo) enters the dynamical 
conductivity o~ij{uj) = a,b,c axis in HTS systems) 
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which generally speaking is a tensor quantity given by 
the following formula 



(Tij(u) = 



d A q 
(2nY 



G(k^q)T J (q,k^q)G(q), 



(62) 



where q = (q, v) and k = (k = 0, uj) and the bare current 
vertex 7^(g, k + q\ k = 0) is related to the Fermi velocity 
vp,i, i.e. Ji(q,k + q; k = 0) = vp^. The vertex function 
Tj(q, k + q) takes into account the renormalization due to 
all sc attering processes responsible for finite conductivity 
jl24| . In the following we study only the in-plane conduc- 
tivity at k = 0. The latter case is realized due to the long 
penetration depth in HTSC cuprates and the skin depth 
in the normal state are very large. In the EPI theory, 
Tj(q, k + q) = Tj(q, iuo ni iuo n + iuj m ) is a solution of an 
approximative integral equation written in the symbolic 
form [68[ 



VeffGGTj 



(63) 



where the effective potential V e ff (due to EPI) is given 
by Veff — I #£ en | 2 ^k, where D K is the phonon 
Green's function. In such a case the Kubo theory predicts 

<T% tra (fj) (i = x,y,z) 



{ dvth{^—)S-\u J ,v) 



f 

Jo 



dv[th(- 



2T 



2T 



(64) 



where S(u, v) = co + S* r (cj + v) - E tr (i/) + i-y l ™ p , 



and 



j tr p is the impurity contribution. In the following we 



omit the tensor index ii in au(u). In the presence of 
several bosonic scattering processes the transport self- 
energy Yi tr (uo) = ReYi tr (uo) + iImYi tr (ou) is given by 



poo 

S tr (w) = - V / di/a? Pi ,Fi(i/)[Jfi(w, i/) + iif 2 (w, i 
, Jo 



(65) 



i^, i/) = + i^) - *(i + i^)], (66) 



K 2 (u,u=^[2cth(^) 



th(^)+th(^)] (67) 



2T 



2T 



Here a^ rl Fi(u) is the transport spectral function which 
measures the strength of the Z-th (bosonic) scattering 
process and ^ is the digamma function. The index I 
enumerates EPI, charge and spin-fluctuation scattering 
processes. Like in the case of EPI, the transport bosonic 
spectral function o% r l F(£l) is given explicitly by 



i 



(v 



x \<l\<l'1 



F(co). 



(68) 



We stress that in the phenomenological SFI theory 



one assumes a 



7V(/i)^ 2 / Im X (k 



p, cj), which, 



as we have repeated several times, can be justified only 
for small g s f, i.e. g s f <C Wb (the band width). 

In case of weak coupling (A < 1), ct(uj) can be written 
in the generalized (extended) Drude form as discussed in 
Section III.B. 
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